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Abstract 



We push further our study of the all-order asymptotic expansion in /3 matrix models with a 
confining, offcritical potential, in the regime where the support of the equilibrium measure is a reunion 
of segments. We first address the case where the filling fractions of those segments are fixed, and show 
the existence of a 1/A^ expansion to all orders. Then, wc study the asymptotic of the sum over filling 
fractions, in order to obtain the full asymptotic expansion for the initial problem in the multi-cut 
regime. We describe the application of our results to study the all-order small dispersion asymptotics 
of solutions of the Toda chain related to the one hermitian matrix model (/3 = 2) as well as orthogonal 
polynomials outside the bulk. 
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1 Introduction 



This paper deals with the ah-order asymptotic expansion for the partition function and muhihnear 
statistics of /3 matrix models. These laws represent a generalization of the joint distribution of the 
N eigenvalues of the Gaussian Unitary Ensemble []\Ieh04] . The convergence of the empirical measure 
of the eigenvalues is well-known (see e.g. [d]\IPS95]), and we are interested in the all-order finite size 
corrections to the moments of this empirical measure. This problem has received a lot of attention in 
the regime when the eigenvalues condensate on a single segment, usually called the one-cut regime. In 
this case, a central limit theorem for linear statistics has been proved by Johansson [Joh98], while a full 
1/iV expansion was derived first for /3 = 2 [APSOl, EM03], then for any /3 > in [BGll]. On the other 
hand, the multi-cut regime remained poorly understood at a rigorous level until recently, except for 
(3 = 2 which is related to integrable systems, and can be treated with the powerful asymptotic analysis 
techniques for Riemann-Hilbert problems, see e.g. [DKM+99b]. Nevertheless, a heuristic derivation 
of the asymptotic expansion for the multi-cut regime was proposed to leading order by Bonnet, David 
and Eynard [BDEOO], and extended to all orders in [Eyn09], in terms of Theta functions and their 
derivatives. It features oscillatory behavior, whose origin lies in the tunneling of eigenvalues between 
the different connected components of the support. These heuristics, initially written for j3 = 2, 
trivially extend to /3 > 0, see e.g. [Borll]. 

Lately, M. Shcherbina has established this asymptotic expansion up to terms of order 1 [Shell, 
Shcl2]. This allows for instance the observation that linear statistics do not always satisfy a central 
limit theorem (this fact was already noticed for /3 = 2 in [Pas06]). In this paper, we go beyond the 
0(1) and put the heuristics of [Eyn09] to all orders on a firm mathematical ground. As a consequence 
for P = 2, we can establish the full asymptotic expansion outside of the bulk for the orthogonal 
polynomials with real-analytic potentials, and the all-order asymptotic expansion of certain solutions 
of the Toda lattice in the continuum limit. The same method would allow to justify rigorously the 
asymptotics of skew-orthogonal polynomials {(3=1 and 4) outside of the bulk, derived heuristically 
in [EynOl]. 

1.1 Definitions 

We consider the probability measure ^ on given by: 

1 ^ 

d/^PJW=^nd^^lB(A.)e-^^(^-) n l^^-^.f- (1-1) 

^Af,,3 i = l lt^i<jt^N 

B is a reunion of closed intervals of Mu{+go}, /3 is a positive number, and Z^'^ is the partition function 
so that (1.1) has total mass 1. This model is usually called the (3 ensemble [Meh04, DE02, ForlO]. We 
introduce the unnormalized empirical measure Afjv of the eigenvalues: 

N 

Mjv = 2'5a.' (1-2) 
and we consider several types of statistics for Mjv. We sometimes denote A = diag(Ai, . . . , Aat). 
Correlators 

We introduce the Stieltjes transform of the n-th order moments of the empirical measure, called 
disconnected correlators: 

Wr,{xi,...,Xn) = fl^'f^ [ / 1 ■ (1.3) 
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They are holomorphic functions of Xi e C\B. For reasons related to concentration of measures, it is 
more convenient to consider the correlators to study large N asymptotics: 



ti=0 



V 

' 1=1 



By construction, the coefficients of their expansion as a Laurent series in the variable Xi — > co give the 
Ti-th order cumulants of A/jv- If / is a set, we introduce the notation xi = (xi)i,=i for a set of variables 
indexed by /. The two type of correlators are related by: 

n s 
s—1 JiO---OJs~I i—1 

If ipn is an analytic function in n variables in a neighborhood of B", the n-linear statistics can be 
deduced as contour integrals of the disconnected correlators: 



^^N!|[ E 'Pn(An,-.-,A,j] = (^^---^^V'n(Cl,.-.,Cn)W^n(Cl,---,Cn). (1.6) 

We remark that the knowledge of the correlators for a smooth family of potentials {Vt)t determines 
the partition function up to an integration constant, since: 

Kernels 

Let c be a n-uple of non zero complex numbers. We introduce the n-kernels: 



K ( \ — ^'^ 



Y\ def^^ {xj - A) 
b=i 



yV;B 

N,0 



(1.8) 



When Cj are integers, the kernels are holomorphic functions of Xj e C\B. When Cj are not integers, the 
kernels are multivalued holomorphic functions of xj in C\B, with monodromics around the connected 
components of B and around co. 

In particular, for (3 = 2, Ki i{x) is the monic A^-th orthogonal polynomial associated to the weight 
lB(a;) ^^^Mx on the real line, and Ar2,(i,-i)(a;, y) is the iV-th Christoffel-Darboux kernel associated 
to those orthogonal polynomials, see Section 2. 

1.2 Equilibrium measure and multi-cut regime 

By standard results of potential theory, see [Joh98] or the textbooks [Dci99, Theorem 6] or [AGZIO, 
Theorem 2.6.1 and Corollary 2.6.3], we have: 

Theorem 1.1 Assume that V : B — > M is a continuous function, and if tod e B, assume that: 

liminf-^^>l. (1.9) 

x^roo 2 In X ^ ^ 
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// V depends on N , assume also that V — > V^^^ in the space of continuous function over B for the 
sup norm. Then, the normalized empirical measure Ln = N^^ Mn converges almost surely and in 
expectation towards the unique probability measure /j,eq := M^'^ on B which minimizes: 

E[^i\ = I dfim^^'HO - II MOMv) In le - v\- (1-10) 

/icq has compact support, denoted S . It is characterized by the existence of a constant C such that: 

VxeB, 2 /d/ieq(01n|a;-e|-^<°^(a;) (1.11) 
Jb 

with equality realized /icq almost surely. Moreover, ifV^'^^ is real-analytic in a neighborhood o/B, the 
support consists of a finite disjoint union of segments: 

S=US,, S, = K,a+], (1.12) 

h=0 

ficq has a density of the form: 

/i=0 

where p'-^ is +1 (resp. —1) if the corresponding edge is soft (resp. hard), and S is analytic in a 
neighborhood ofS. 

The goal of this article is to establish an all-order expansion of the partition function, the correlators 
and the kernels, in all such situations. 

1.3 Assumptions 

We will refer throughout the text to the following set of assumptions. 
Hypothesis 1.1 

• (Regularity) V : B — > M is continuous, and ifV depends on N, it has a limit V^^^ in the space 
of continuous functions over [6_, &+] for the sup norm. 

• (Confinement) If too e B, liniinfa;^^,3o ^^^^^i > 1- 

• (g + l-cut regime) The support of is of the form S = U^=o where Sh = [ck;j^,Q;,^] with 

• (Control of large deviations) The effective potential [/^'^(x) = V{x) — 2 Jin jx — V\dlJ^^{C} 
achieves its minimum value on S only. 

• (Offcriticality) In the equilibrium measure (1.13), S{x) > m S. 

At some point, we shall need to add a stronger assumption concerning off-criticality: 
Hypothesis 1.2 The same as Hypothesis 1.3, and: 

• (Strong off-criticality) For any soft edge a*, we have 5"(a^) ^ 0. 
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We will also require regularity of the potential: 
Hypothesis 1.3 

• (Analyticity) V extends as a holomorphic function in some open neighborhood U o/S. 

• (1/N expansion of the potential) There exists a sequence {V^''^)k^o of holomorphic functions in 
U and constants (w^'^^)fe>o such that, for any K ^ 0, 



sup 



K 

^(0 - 2 v<^5 7v-(^+i). (1.14) 

fc=0 

In Section 6, we shall weaken Hypothesis 1.3 by allowing complex perturbations of order 1/A^ and 
harmonic functions instead of analytic functions; 

Hypothesis 1.4 y : B ^ C can be decomposed as V = Vi + V2 where: 

• For j = 1,2, Vj extends to a holomorphic function in some neighborhood U o/B. There exists a 
sequence of holomorphic functions (vj'^^)fe^o o,iT'd constants (w]'^^)fc>o so that, for any K > 0.' 



sup 



K 



k=0 



^AO - 2 N-''vj''\0 < 7V-(^+i). (1.15) 



• l/{0} = v{ ' + V2* ' IS real-valued on B. 

The topology for which we study the large N expansion of the correlators is described in § 5.1, and 
amounts to controlling the (moments of order m) x C™ uniformly in m for some constant C > 0. We 
now describe our strategy and announce our results. 

1.4 Main result with fixed filling fractions 

Before coming to the multi-cut regime, we analyze a different model where the number of A's in a 
small enlargment of is fixed. Let A = Uh=o^'» where = [a/^,a/|^] are pairwise disjoint segments 
such that < aj^ < ctt < ■ introduce the set: 

f ^ 1 

£g = [ee]0,l[', <l}- {1-16) 

/i=i 

If e e Eg, we denote eo = 1 ^ l^h=i^h' we let TV = ([iVeoJ, [iVeiJ, . . . , [A^egJ), and consider the 

l/i=0 

g ^ Nh 



probability measure on nfi=o ^h'^ ■ 



-1 y ^vfa 

'^N,e,l3 h=0 i=l lt^i<jt^N 



X 



n n iAM-A.',,f. (1.17) 

The empirical measure Mjv,e and the correlators Wn,eixi, . . . ,Xn) for this model are defined as in 
§ 1.1. We call eh the filling fraction of Ah. It follows from the definitions that: 

WnA^,X2,...,Xn)=Sn,lNeh. (1.18) 



A. 2i7r 



We will refer to (1.1) as the initial model, and to (1.17) as the model with fixed filling fractions. 
Standard results from potential theory imply: 
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Theorem 1.2 Assume V regular and confining on A. Then, the normalized empirical measure 
N^^ Mpf,^ converges almost surely and in expectation towards the unique probability measure /icq.e 
on A which minimizes: 

E[t^] = II (^(W°>(a + W°>(77))-ln|e-7y|)d/i(e)dM^). (1.19) 

among probability measures with partial masses /i[A^] = Ch- They are characterized by the existence 
of constants C^.h such that: 

\fxeAh, 2 f d^i,^.,{Olli\x-^\-V^°Hx)^C,,h, (1.20) 
Jb 

with equality realized /icq.c almost surely. fJ-cq.e can be decomposed as a sum of positive measures /icq.e,/i 
having compact support in A^, denoted Sgji. Moreover, if V^^^ is real-analytic in a neighborhood of 
A, Se^/i consists of a finite reunion of segments. 

/icq appearing in Theorem 1.1 coincides with /icq.e, for the optimal value e* = (/icq[Aft])i<;ft^g, and in 
this case Se^,,h is actually the segment [a^,Q;,'|^]. The key point is that, for e close enough to e^,, the 
support Se.h remains connected, and the model with fixed filling fi:action enjoys a 1/N expansion. 

Theorem 1.3 IfV satisfies Hypotheses 1.1 and I.4 on A, there exists t > such that, uniformly for 
e e f g such that |e — e*| < t, we have an expansion for the correlators: 

M^„,,(a;i,...,x„) = ^ 7V-^-M^W(a.i,...,a;„) + 0(7V-=^). (1.21) 

Up to a fixed 0{N^^) and for a fixed n, (1.21) holds uniformly for xi, . . . ,Xn in compact regions of 
C\A. Besides, if the strong off-criticality of Hypothesis 1.2 is satisfied, wii^} are smooth functions of 
e close enough to e». 

We prove this theorem, independently of the nature soft/hard of the edges, in Section 5 with real- 
analytic potential (i.e. Hypothesis 1.3 instead of 1.4). The result is extended to harmonic potentials 
(i.e. Hypothesis 1.4) in Section 6.2. Actually, we provide in Proposition 5.5 an explicit control of the 
errors in terms of the distance of xi, . . . , to A, and its proof makes clear that the expansion of the 
correlators is not expected to be uniform for xi, . . . ,Xn chosen in a compact of C\A independently of 
n and K. 

We then compute in Section 7.1 the expansion of the partition function thanks to the expansion 
of Wi^e, by a two-step interpolation preserving Hypotheses 1.2-1.4 between our potential V and a 
reference situation where the partition function is exactly computable for finite TV, in terms of Selberg 
integrals. 

Theorem 1.4 IfV satisfies Hypotheses 1.2 and I.4 on A, there exists t > such that, uniformly for 
e e £g such that |e — e*| <t, we have: 



y ZnIIp = exp ( 2 iV-^ 45 + 0{N'^n) , (1-22) 



with e = y,l-r,s„- „+, where: 

3 + /3/2 + 2//3 /3/2 + 2//? -1 + 2//3 + /3/2 
e++ = ' " " 6 ' " 4 ' ^ 

and we recall p* = 1 for a soft edge and p* = —1 for a hard edge. Besides, F^^^ is a smooth function 
of e close enough to e^, , and at the value € = the derivative of ^ vanishes and its Hessian is 
negative definite. 
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Up to a given 0{N~^), all expansions are uniform with respect to parameters of the potential and 
of e chosen in a compact set so that the assumptions hold. The power of N in prefactor is universal 
in the sense that it only depends on the nature of the edges, and its value can be extracted from the 
large N expansion of Selberg type integrals. Theorems 1.3-1.4 are the generalizations to the fixed 
filling fraction model of our earlier results about existence of the 1/N expansion in the one-cut regime 
[BGll] (see also [Joh98, APSOl, EM03, GMS07, KSIO] for previous results concerning the one-cut 
regime in /3 = 2 or general /? ensembles). 

1.5 Main results in the multi-cut regime 

Let us come back to the initial model (1.1), and take A = IJ^^q A/i c B a small enlargement of 

V:B 

N,I3 



the support S as in the previous paragraph. It is well-known that the partition function zYr'^ can 



be replaced by z'^''^ up to exponentially small corrections when N is large (see [] for results in this 
direction, and we give a proof for completeness in § 3.1 below). The latter can be decomposed as a 
sum over all possible ways of sharing the A's between the segments A^, namely: 

where we have denoted Nq = N — Xi/^=i -^h the number of A's put in the segment Aq. So, we can use 
our results for the model with fixed filling fractions to analyze the asymptotic behavior of each term 
in the sum, and then find the asymptotic expansion of the sum taking into account the interference 
of all contributions. 

In order to state the result, we need to introduce the Siegel Theta function with characteristics 
/X, 1/ 6 C^. li T he a g X g matrix of complex numbers such that Imx > 0, it is the entire function of 
V e defined by the converging series: 



1? 



{v\t) = ^ exp (m{rn + ^) ■ T ■ (m + fi) + 2m{v + u) ■ (m + fjb)] . (1-25) 



Among its essential properties, we mention: 

• for any characteristics /.t, f , it satisfies the diffusion-like equation Aindrf^ h''^ ~ ^'"h^'u^i'&- 

• it is a quasiperiodic function on the lattice ©t(Z3): for any mo, no e 

^"\v\t). (1.26) 



{v + mo + T ■ tiqIt) = exp (2i7rmo • /x — 2i7rno ■ {v + v) — ivrno • r • no) i9 



• it has a nice transformation law under r {At + B){Ct + D) ^ where A, B, C, D are the 

g X g blocks of a 23 x 2g symplectic matrix [Mum84]. 

• when T is the matrix of periods of a genus g Riemann surface, it satisfies the Fay identity [Fay 70] . 

We define the operator V^, acting on the variable v of this function. For instance, the diffusion 
equation takes the form 4i7r(?T-i? = V®^z9. 

Theorem 1.5 Assume Hypotheses 1.2 and 1.4- Let e* = (/icq[Sh])i^hs;g- Given the coefficients of 
the expansion in the fixed filling fraction model from Theorem I.4, we denote (-pf^"^)^^-* their tensor of 
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£-th order derivatives with respect to e, evaluated at e* . Then, the partition function has an asymptotic 
expansion of the form: 



yV■.^ 







In this expression, if X is a vector with g components, T^^}\X\ = 1, and for 1: 



r 



>9 

mi,.. .,mr^ — 2 
Y,l_iei+m,=k 



where ■ denotes the contraction of tensors. We have also introduced: 



(Ft;')' 

2i7r ' 



2m 



Being more explicit but less compact, we may rewrite: 







mi ,. . . ,m,.^— 2 



2 (m-iVe.)®^^ 



(1.27) 



(1.28) 



(1.29) 



(1.30) 



For /3 = 2, this result has been derived heuristically to leading order in [BDEOO], and to all orders 
in [Eyn09], and the arguments there can be extended straightforwardly to all values of /3, see e.g. 
[Borll]. Our work justifies their heuristic argument. We exploit the Schwinger-Dyson equations for 
the /3 ensemble with fixed filling fractions taking advantage of a rough control on the large N behavior 
of the correlators. The result of Theorem 1.5 has been derived up to o(l) by Shchcrbina [Slicl2] for real- 
analytic potentials, with different techniques, based on the representation of Oiss/K/i'sSg — ^h'.j \^ , 
which is the exponential of a quadratic statistic, as expectation value of a linear statistics coupled 
to a Brownian motion. The rough a priori controls on the correlators do not allow at present the 
description of the o(l) by such methods. The results in [Shcl2] were also written in a different form: 
F^"^ was identified with a combination of Fredholm determinants (see also the physics paper [WZ06]), 
whereas this representation does not come naturally in our approach) . Also, the step of the analysis 
of Section 8 consisting in replacing the sum over nonnegative integers such that Nq + . . . + Ng = N in 
(1.24), by a sum over AT e Z^, thus reconstructing the theta function, was not performed in [Shcl2] 

Let us make a few remarks. The 2i7r appears because we used the standard definition of the Siegel 
theta function, and should not hide the fact that all terms in (1.30) are real-valued. Here, the matrix: 



Hessian(Fj;g^^) 
2i7r 



(1.31) 



involved in the theta function has purely imaginary entries, and Imr*^^ is definite positive according 
to Theorem 1.4, hence the theta function in the right-hand side makes sense. Notice also that for 
it is Z^-periodic in its characteristics fi, hence we can replace — iVe» by —Ne^ + [A^e,J, and this is 
responsible for modulations of frequency 0{1/N) in the asymptotic expansion, and thus breakdown of 
the 1/N expansion. Still, "subsequential" asymptotic expansions in 1/A'' may occur. For instance, in 
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a symmetric two cuts {g = 1) model, we have e* = 1/2 and thus the right-hand side is an asymptotic 
expansion in powers of 1/N depending on the parity of N. 
Let us give the two first orders of (1.30): 



and: 



T*^j[X] = I iF}-'Y • + 1 (f{-i>)" . + (Fi5)' . X, (1.32) 



Tl'^iX] = ^[(Fi/)'1®'-X®6+ j^[(Fi/)'''®(Fi:/^)''].X®^ 



+ (2 [K3)f' + 2 ■ + • (1-33) 

If the potential V is independent of /3, we observe that e* does not depend on /3, and it is well-known 
[CE06] that the coefficients in the expansion (1-22) have a simple dependence in /3: 

[fc/2j + l ^ 1-G 

In particular, those coefhcients vanish for odd k. The first few ones are: 



S f) (l-|)'""%i««-l. (1.34) 



G=0 



= Ft;' = - 1) ^["^1], = ^[1-"] + (f + I - 2) (1.35) 

and have been first identified in [WZ06] . In particular, for the argument of the theta function: 

"-^=u-^j^i^' ^-''=2^i^- (^-^^^ 

Similarly for the correlators in the fixed filling fraction model, the dependence in /3 takes the form: 

l(fe-»i+2)/2j . . l-G-n 0,/c + 2-2G-n 

Wi'h^i,-.-,^n)= E (f) (l-|) We'^-^^-"](.i,...,cr„). (1.37) 

All coefficients fJ*^'^! and functions Wn^e^\xi, . . . ,a;„) can be computed with the /? deformation of 
the topological recursion formulated by Chekhov and Eynard [CE06], applied to the spectral curve 
determined by the equilibrium measure /icq,e and W2 which encodes the covariance of linear statistics 
at leading order in the model with fixed filling fractions. We stress now a point of this theory relevant in 
the present case. When is a polynomial and e is close enough to e«, the density of the equilibrium 
measure can be analytically continued to a hyperelliptic curve of genus (7, denoted Cg and called 
spectral curve. Its equation is: 

y'=Yl{x-a-^ri{x-atjy^. (1.38) 

Let be the cycle in Ce surrounding Agji = The family ^ = (^^)i^;i^j, can be completed 

by a family of cycles B so that {A., B) is a symplectic basis of homology of C^. The correlators wk%'^^ 
are meromorphic functions on C^, computed recursively by a residue formula on Cg. In particular, the 
analytic continuation of 

0/ \ ia;[o.o]/ n , , 2 da;idx2 , 
u}2{xi,X2) = W2 e (a;i,a;2)da;ida;2 + -5 -, (1.39) 

P [xi — X2) 
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is the unique 2-forni on C^, which has vanishing A. periods, and has for only singularity a double 
pole with leading coefficient ^ and without residue at coinciding points. Then, it is a property of 
the topological recursion that the derivatives of F^'^'^^ can be computed as ,B-cycle integrals of the 
correlators: for any {G,K) ^ (0,0), (0,1), 

(j-[G,K])W = i dx^...(f) Ax,wff\xi,...,x,), (1.40) 

is JB 

and for any (n, G, K): 

(W,P'^])(^)(a;i, ...,Xn) = (f dx„+i ■■■(() dx^+e W,!,+fJ(a:i, . . . , Xn+e). (1.41) 

Jb Jb 

In particular: 

{Wff)'{x)dx = 2{ttv^{x) (1.42) 

where -uj is the basis of holomorphic 1-forms on dual to A., i.e. characterized by <^^^ Wh' = 5h.h'- 
This formula at e = e» can be used to compute the functions Tg'^^[X] appearing in (1.28). The 
derivation wrt e is not a natural operation in the initial model when N is finite, since Neh are forced 
to be integers in (1.17). We rather show that the coefficients of expansion themselves are smooth 
functions of e, and thus de makes sense. 

For /3 = 2, we remark from (1.34) that the coefficients Pl'^j^^^^ all vanishes, so that we retrieve 
the celebrated 1/iV^ expansion in the one-cut regime or in the fixed filling fraction model. This is 
in general not true anymore in the multi-cut regime. For instance, we have a term of order 1/A^ 
involving: 

tI%,[X] = \ (F<-^>)"' . X®3 + (FiS)' . X. (1.43) 

In a two-cut regime {g = 1), a sufficient condition for all terms of order 7V~*^^'^+'^^ to vanish is that 
= 1/2 and Z^''^ = Z^'^_^, i.e. the potential has two symmetric wells. In this case, we have an 
expansion in powers of 1/7V^ for the partition function, whose coefficients depends on the parity of 
N. In general, we also observe that t'*.^=2 = 0, i.e. Thetanullwerten appear in the expansion. 

1.6 Asymptotic expansion of kernels and correlators 

Once the result on large N expansion of the partition function is obtained, we can easily infer the 
asymptotic expansion of the correlators and the kernels by perturbing the potential by terms of order 
1/iV, maybe complex- valued, as allowed by Hypothesis 1.4. 

1.6.1 Leading behavior of the correlators 

Although we could write down the expansion for the correlators as a corollary of Theorem 1.5, we 
bound ourselves to point out their leading behavior. Whereas Wn behaves as 0{N'^~^) in the one-cut 
regime or in the model with fixed filling fractions, Wn for n 5= 3 does not decay when N is large in a 
{g + l)-cut regime with g ^ I. More precisely: 

Theorem 1.6 Assume Hypothesis 1.2, I.4 and number of cuts (g + 1) ^ 2. We have, for uniform 
convergence when xi, . . . , x„ belongs to any compact of (C\A)"' : 

{v.,f,\T.,^), (1.44) 
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W2ixi,X2) 



N~ 



W^°}{X1,X2) + 



dxi 



dx-2 



• V®2 In I? 







and for any n > 3; 







v^.pW^.p)- (1-45) 



Integrating this result over ^-cycles provide the leading order behavior of n-th order moments of the 
filling fractions TV. We will also describe in Section 8.2 the fluctuations of the filling fractions: we 
find that they converge to a discrete Gaussian random variable. 

1.6.2 Kernels 

We explain in § 0.4 that the following result concerning the kernel - defined in (1.8) - is a consequence 
of Theorem 1.3: 

Corollary 1.7 Assume Hypothesis 1.2 and 1.4. There exists t > such that, for any e e Eg such 
that |e — e*| < t, the n-kernels in the model with fixed filling fractions have an asymptotic expansion 
of the form: 

if„,c,.(a:i,...,x„) =exp[ 2 N-'^{j]l,/^tMr\':h)+0{N-^)l (1.46) 

k^-l n=l 

where Cx^c is the linear form : 

>Cx,c = E / . (1.47) 

Up to a given 0{N~^), this expansion is uniform for xi, . . . ,Xn in any compact o/C\A. 
Hereafter, if 7 is a smooth path in C\Se, we set £j = J^, and £®" is given by: 

£®"[M/W]= [dx^--- f dx„M^W(a;i,...,x„). 



A priori, the integrals in the right-hand side of (1.4G) depend on the homology class in C\A of paths 
CO — > Xi. A basis of homology cycles in C\A is given by ^ = {Ah)o^h^g- We also denote for 
convenience e = (eh)os;/is;g- We deduce from (1.18) that: 

^Wi''}{^,X2, . . . ,Xn) = Sn,lSk,-ie. (1.48) 

Therefore, the only multivaluedness of the right-hand side comes from the first term N J d£,Wi ^^\^), 
and given (1.48) and observing that Ne^ are integers, we see that it exactly reproduces the mon- 
odromies of the kernels depending on Cj. 

We now come to the multi-cut regime of the initial model. If X is a vector with g components, 
and £ is a linear form on the space of holomorphic functions on C\S£, let us define: 

= tv^ E ^^ ^"^mj^fm ^ . ^^(ZL...)^ (1.49) 

mi ,. . . ,771^.^—2 

ni ,.. .,71^,^0 

where we took as convention W^lo g = Fe^K Then, as a consequence of Theorem 1.5: 
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Corollary 1.8 Assume Hypothesis 1.2 and 1.4- With the notations of Corollary 1.7, the n-kernels 
have an asymptotic expansion^: 

if„,c(x) = if„,e,*(x)(l + 0(7V-*')) (1.50) 
























where Cx,c = ''j l^J ^'^'^ ba,sis of holomorphic 1-forms. 

A diagrammatic representation for the terms of such expansion was proposed in [BE12, Appendix A]. 

2 Application to orthogonal polynomials and integrable sys- 
tems 

Since orthogonal polynomials are related to the 1-hcrmitian matrix model (i.e. (3 = 2), our results can 
be used to establish the all-order asymptotics of orthogonal polynomials outside the bulk (Theorem 2.2 
below). We will illustrate it for orthogonal polynomials with respect to an analytic weight defined on 
the whole real line, but it could be applied equally well to orthogonal polynomials with respect to an 
analytic weight on a finite union of segments of the real axis. 

The leading order asymptotic of orthogonal polynomials is well-known since the work of Deift et 
al. [DKM+97, DKM+99b, DKM+99a], using the asymptotic analysis of Riemann-Hilbert problem 
which was pioneered in [DZ95]. In principle, it is possible to push the Riemann-Hilbert analysis 
beyond leading order, but this approach being very cumbersome, it has not been performed yet to 
our knowledge. Notwithstanding, the all-order expansion has a nice structure, and was heuristically 
derived by Eynard [Eyn06] based on the general works [BDEOO, Eyn09]. In this article, we provide a 
proof of those heuristics. 

Unlike the Riemann-Hilbert technique which becomes cumbersome to study the asymptotics of 
skew-orthogonal polynomials (i.e. /3 = 1 and 4) and thus has not been performed up to now, our 
method could be applied without difficulty to those values of /3, and would allow to justify the heuristics 
of Eynard [EynOl] formulated for the leading order, and describe all subleading orders. In other words, 
it provides a purely probabilistic approach to address asymptotic problems in integrable systems. It 
also suggest that the appearance of theta functions is not intrinsically related to integrability. In 
particular, we see in Theorem 2.2 that for /3 = 2, the theta function appearing in the leading order is 
associated to the matrix of periods of the hyperelliptic curve C^^ defined by the equilibrium measure. 
Actually the theta function is just the basic block to construct analytic functions on this curve, and 
this is the reason why it pops up in the Riemann-Hilbert analysis. However, for /? 2, the theta 
function comes is associated to (/3/2) times the matrix of periods of C^^, which might be or not the 
matrix of period of a curve, and anyway is not that of C^^^ . So, the monodromy problem solved by this 
theta function is not directly related to the equilibrium measure, which makes for instance for (3=1 
or 4 its construction via Riemann-Hilbert techniques a priori more involved. 

Contrarily to Riemann-Hilbert techniques however, we are not yet in position within our method 
to consider the asymptotic in the bulk, at the edges, or the double-scaling limit for varying weights 
close to a critical point, or the case of complex- values weights which has been studied in [BM09]. We 
hope those technical restrictions to be removable in a near future. 

-"^We warn the reader that ' denotes a derivative with respect to filling fractions, not with respect to variables of the 
correlators. 
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2.1 Setting 

We first review the standard relations between orthogonal polynomials on the real hne, random ma- 
trices and integrable systems, see e.g. [CG12, Section 5]. In this section, /3 = 2 and we omit to 
precise it in the notations. Let T4(A) = V{X) +X;fe=i ^k^''- Let {Pn,N{x))n^o be the monic orthogonal 
polynomials associated to the weight dw{x) = da; e"^^*'-^-' on B = M. We choose V and restrict in 
consequence so that the weight decreases quickly at +go. If we denote hn^N the L^{dw) norm of 
Pn.N, the polynomials Pn.N = Pn,N / \/hn,N are orthonormal. They satisfy a three-term recurrence 
relation: 



xPn^N{.x) = ^/K^Pn+i{x) + (3nM Pnix) + y'/l„_i, ArP„_ i (x) . (2.1) 

The recurrence coefficients are solutions of a Toda chain: if we set 

Un.N = ln/l„,Ar, Vn.N = —Pn.N, (2.2) 

we have: 

3tiUn,N = Vn,N - Vn^l,N, Sf^Vn^N = C^-'+i'" - e""'", (2.3) 

and the coefficients tk generate higher Toda flows. The recurrence coefficients also satisfy the string 
equations: 

ri 

^/^h^[V'{QN)]n,n^l = ^, [V'{QN)]n,n = 0, (2.4) 

where Qn is the semi-infinite matrix: 



Q 



N 



V J 



(2.5) 



Eqns 2.4 determines in terms of V the initial condition for the system (2.3). The partition function 
T(t) = Z'^'^ is the Tau function associated to the solution (un,jv(t), fn,w(t))„^i of (2.3). The 
partition function itself can be computed as [Mch04, PSll]: 

Zn^ = ^! n ''i-^- (2.6) 

We insist on the dependence on and V by writing hj^N = hj{NV). Therefore, the norms can be 
retrieved as: 



n" h(NV) 1 r^NV/(n+l)-R + 

UUh.iNV) n + 1 n + 1 zV''^ ' ^' 

The regime where n, — > oo but s = njN remains fixed and positive correspond to the small 
dispersion regime in the Toda chain, where 1/n plays the role of the dispersion parameter. 

2.2 Small dispersion asymptotics of /i„ 7v 

When Vto/so satisfies Hypotheses 1.2 and 1.3 for a given set of times (so,to), Vtl s satisfies the same 
assumptions at least for (s,t) in some neighborhood U of (sojto), and Theorem 1.5 determines the 
asymptotic expansion of T/v(t) = Z^"^ up to 0(7V~^). Besides, we can apply Theorem 1.5 to study 
the ratio in the right-hand side of (2.7) when n ^ oo. 
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Theorem 2.1 In the regime n, N ^ co, s = n/N > fixed, and Hypotheses 1.2 and 1.3 are satisfied 
with soft edges, we have the asymptotic expansion: 



Un,N 



+ ln 



-{n + 1) e* 








(o|n] 



2-2G-m<0 



+ lll 



-In 



1 + 



-(n + 1) e« 




1 + 



V 







-{n 


+ l)e/ 







—71 e, 



(o|n) 


-& 





(o|n) 



Here, are the filling fractions of fjXq'^ and Cy^/g is the linear form defined by: 



S 2i7r s 



(2.8) 



(2.9) 



We have not performed the expansion of l/(n + 1) in powers of 1/n to make the structure more 
transparent. We recaU that all the quantities Wm,* can be computed from the equilibrium measure 
associated to the potential Vt, so making those asymptotic explicit just requires to solve the scalar 
Riemann-Hilbert problem for /i^^'. Notice that the number g + 1 of cuts a priori depends on (so,to), 
and we do not address the issue of transitions between regimes with different number of cuts (because 
we cannot relax at present our off-criticality assumption), which are expected to be universal [Dub08]. 
We also collect here in one place and for /? = 2 other notations appearing throughout the text: 



^IG] ^ ^[G] ^ p{2G-2} y^;[G] ^ ^{2G- 



2+n} 



(jfV 

2i7r ' 



and 



fW[/:;X] = 



-I 



[G.b(£0 



Gi,...,Gr>0 
ii + 2Gi-2>0 



(X) 

1=1 



(JQ 

£■! 



(2.10) 



(2.11) 



-I 



fl,...,fr>l 

Gi,...,G^>0 

ni ,. . .jTlr^O 

£i+2Gi-2+ni>a 
Sr_i(^. + 2G.-2+n.) = fc 



/:®"'[(wif,;i)(^')] 



.1 /.I 



(2.12) 
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2.3 Asymptotic expansion of orthogonal polynomials away from the bulk 

The orthogonal polynomials can be computed thanks to Heine formula [Szc39]: 

n 

P^{x) = /ij:'/^*[n(^ - ^^)] = ^i^i(^)- (2.13) 

1=1 

Hence, as a corollary of Theorem 1.8: 

Theorem 2.2 In the regime n, N ^ co, s = n/N > fixed, and Hypotheses 1.2 and 1.3 are satisfied, 
for X e C\S, we have the asymptotic expansion: 

P„(x) = exp ( ^ ^ „2-2G-„. ^r"N3] A (1 ^ (2 14) 























(0 


-0 



where Lx = j^,. Up to a given 0{N this expansion is uniform for x in any compact o/C\S. 

We remark that >Ca;[-cc7] is the Abel map evaluated between the points x and oo. 

As such, the results presented in this article do not allow the study of the asymptotic expansion 
of orthogonal polynomials in the bulk, i.e. for x e S. Indeed, this requires to perturb the potential 
V{X) by a term — ^ ln(A — x) having a singularity at x e S, a case going beyond our Hypothesis 1.4. 
Similarly, we cannot address at present the regime of transitions between a g cut regime and and 
g'-cut regime with g ¥= g' , because offcriticality was a key assumption in our derivation. Although it 
is the most interesting in regard of universality, the question of deriving uniform asymptotics, even at 
the leading order, valid for the crossover around a critical point is still open from the point of view of 
our methods. 



3 Large deviations and concentration of measure 
3.1 Restriction to a vicinity of the support 

Our first step is to show that the interval of integration in (1.1) can be restricted to a vicinity of 
the support of the equilibrium measure, up to exponentially small corrections when N is large. The 
proofs are very similar to the one-cut case [BGll], and we remind briefly their idea in § 3.2. Let V 
be a regular and confining potential, and /i^'^ the equilibrium measure determined by Theorem 1.1 
or Theorem 1.2. We denote S its (compact) support. We define the effective potential by: 

U'''^{x) = V{x)-2 f d^,Yf{Oln\x-^\, C/^^^(a;) = f/^^^(:r)-inf[/^^^(e) (3.1) 

when a; e B, and +oo otherwise. 

Lemma 3.1 IfV is regular and confining, we have large deviation estimates: for any F c B\S closed 
and c B\S open, 

limsup^lnAipppz A, e F] < inf ;7^-B(x), (3.2) 

liminf 4lnM«'«[3i e 0] 5= inf C7^'^(a;). (3.3) 

N—^oo N 2 xeO 
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We say that V satisfies a control of large deviations on B if C/^'^ is positive on B\S. Note that U^'^ 
vanishes at the boundary of S. According to Lemma 3.1, such a property implies that large deviations 
outside S are exponentially small when N is large. 

Corollary 3.2 Let V be regular, confining, satisfying a control of large deviations on B, and assume 
dB nS = 0. Let A c B 6e a finite union of segments such that S c A. There exists r](A) > so that: 

^]^;' = 4;"(l + 0(e-^''^'^))), (3.4) 

and for any n ^ 1, there exists a universal constant 7„ > so that, for any Xi, . . . , a;„ e (C\B)".' 

= -W»?(A) 

It is useful to have a local version of this result: 



W^-^^'ix,, . . . ,Xn) - W^'-^ixi, . . . ,Xn)\ < Jn \, ■ (3.5) 



Corollary 3.3 Let V be regular, confining, satisfying a control of large deviations on B, and assume 
dB n S = 0. Let A c B 6e a finite union of segments such that S c: A. If is the left (resp. right) 
edge of a connected component of A, let us define A^ = A u [a, ao]. For any e > small enough, there 
exists rj^ > so that, for N large enough and any a ejao — e, oq + e[, we have: 



v■.^a 



da In ZJ 



^6-^"% (3.6) 



and, for N large enough and any n 5= 1 and Xi, . . . , x^ e {C\Aa) : 

da' W^'^^ix^, . . . ,X„) < r (3.7) 

From now on, even though we want initially to study the model on B^, we are going first to study 
the model on A^, where A is small (but fixed) enlargement of S as allowed above, in particular we 
choose A bounded. 

Proposition 3.4 For any fixed e e £g, the same results holds for the partition function and the 
correlators in the fixed filling fraction model. 

3.2 Sketch of the proof of Lemma 3.1 

We only sketch the proof, since it is similar to [BGll]. 

Recall that L^ = -^^^ Xlili '^Ai denotes the normalized empirical measure, either in the initial 
model, or in the fixed filling fraction model. We represent: 

MpJ[3z A.6F] = 7V^^ (3.8) 

where, for any measurable set X: 



V B ivv .g 



dC exp -i - ^ V{0 + {N- 1)13 I diAr_i(A) In |e - A| 



(3.9) 



We first prove a lower bound for T]^'^(X) assuming X contains at least an open interval, of size larger 
than some e > 0. Let ki{V) be the Lipschitz constant for V on B, and: 

X^ = {a;6B, inf la; - CI > e/2} (3.10) 

CeB\X 
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Using twice Jensen's inequality, we get 



NV .R 

Ivr — r 



TV:B /\/\ N-1 



x+e/i 
x—s/i 



2:+£/4 



X — 6/i 



d^exp (iV-l)/3 / dLjv_i(A)ln|^-A| 



> - sup e 



•/ B 



where we have set: 



Hx.eW 



x+e/i 



ln|f-A| 



(3.11) 



(3.12) 



Ix-e/i 

For any fixed e > 0, H^.e is bounded continuous on any compact, so we have by Theorem 1.1 in the 
initial model (or Theorem 1.2 in the fixed filling fraction model): 



T];;;^ (X) 5* I sup e-^ exp [{N - l)/3 d^^f (A) i?,,,(A) + iVi?(£, TV)} (3.13) 



with limjv^x R{£,N) = 0. Letting N ^ co, we deduce: 



(3.14) 



Interchanging the integration over ^ and x, observing that £, ^ J d/iJ^'^(A) In |^ — A| is smooth and 
then letting e ^ we conclude 



liminf — lnT]C'«(X) > MU^'^(x) 

N^OD N ^^^^ ' 2 X€X ^ ' 



(3.15) 



where we have recognized the effective potential of (3.1). To prove the upper bound, we observe that 
for any M > 0, 

TTPJW /^S;^ / deexp{-^V(e)-(iV-l)/? I dL^_i(A)lnmax(|e-AUf-i)} 

(3.16) 

As A ^ lnmin(|^ — A|,M~^) is bounded continuous on compacts, we can use Theorem 1.1 in the 
initial model (or Theorem 1.2 in the fixed filling fraction model) to deduce that for any e > 



T];;^(X) ^ / dC exp 



with 



(3.17) 
(3.18) 



1 JVV .g 

limsupi?(£, TV) = limsup — ln^j!Jl\'o(d(LAr_i,/iJl'^) > e) < 0. 



Moreover, ^ ^(^) ^ / d/i^'^(A) lnniax(|^ — A|,Af ^) is bounded continuous so that a standard 
Laplace method yields 



liniinf 1 In T];;^(X) =S - mf { ^ (v{0 In - A| v Af-^ 



(3.19) 



Finally, we conclude by monotone convergence theorem which implies that / d/iJ^'^(A) lnmax(|^ — 
A|, M~^) increases as M goes to infinity towards / d^^'^(A) In |f — A|. 
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3.3 Concentration of measure and consequences 

We will need rough a priori bounds on the correlators, which can be derived by purely probabilistic 
methods. This type of result first appeared in the work of [dMPS95, Joh98] and more recently 
[KSIO, MMS12]. Given their importance, we find useful to prove independently the bound we need 
by elementary means. 

Hereafter, we will say that a function / : M ^ C is &-Lipschitz if 

Kb(f) = sup — |T — < CO. (3.20) 

Our final goal is to control {Lj^ — ficq)[^] for a class of functions ip which is large enough, in particular 
contains analytic functions on a neighborhood of the interval of integration A. This problem can be 
settled by controlling the "distance" between and /icq for an appropriate notion of distance. We 
introduce the pseudo-distance between probability measures: 

S)(m, u) = - jj d[^Ji- u]{x)A[ti - y]{v) In \x - y\ (3.21) 

which can be represented in terms of Fourier transform of the measures by: 

^{ii,u)= r p-\{^-V){sf (3.22) 



Since Ljq has atoms, its pseudo-distance to another measure is in general infinite. There are several 
methods to circumvent this issue, and one of them, that we borrow from [MMS12], is to define a reg- 
ularized measure (see the beginning of § 3.4.1 below) from L^- Then, the result of concentration, 
takes the form: 

Lemma 3.5 Let V he regular, , confining, satisfying a control of large deviations on A. There 
exists C > so that, for t small enough and N large enough: 

M]^;^[2)[ZiV,Mer] e<^^l"^-^^*^ (3.23) 

We prove it in § 3.4.1 below. The assumption V of class ensures that the effective potential (3.1) 
defined from the equilibrium measure is a ^-Lipschitz function (and even Lipschitz if all edges are 
soft) on the compact set A, as one can observe on (A. 4) given in Appendix A. 
This lemma allows a priori control of expectation values of test functions: 

Corollary 3.6 Let V be regular, C'^ , confining, satisfying a control of large deviations on A. Let 
b> 0, and assume ip : R ^ C is a b- Lipschitz function with constant Kij{(f), and such that: 

|^|i/2 := ( I ds|s||^(s)|2)'^' <(». (3.24) 

Then, there exists C3 > such that, for t small enough and N large enough: 

j^A[LN A^r^^Kx) ^{x)\ ^ (^f^ + t I^K/,] ^ e'^3^^'"~-^^*^ (3.25) 

As a special case, we can obtain a rough a priori control on the correlators: 

Corollary 3.7 Let V he regular, , confining and satisfying a control of large deviations on A. Let 
D' > 0, and: 
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There exists a constant 7i(A, _D') > so that, for N large enough: 

\Wiix)-Nwt'\x)\^^iiA,D')wNf{d{x,A)). (3.27) 
Similarly, for any n > 2, there exists constants 7^i(A,_D') > so that, for N large enough: 

n 

\Wr,{xi,...,Xn)\^ln{A.D')wlY\f{d{x,,k)). (3.28) 



In the {g + l)-cut regime with 5 > 1, we denote {Shjo^^h^^g the connected components of the support 
of , and we take A = U^=o where kh = [a,~, a,^] c B are pairwise disjoint bounded segments 
such that S?i c A^. For any configuration A e A^, we denote the number of Aj's in A/j, and 
AT = {Nii)i^hiig- The foUowing result gives an estimate for large deviations of N away from A^e, in 
the large N limit. 

Corollary 3.8 Let A be as above, and V be , confining, satisfying a control of large deviations on 
A, and leading to a (g + I) -cut regime. There exists a positive constant C such that, for N large 
enough and uniformly in t: 

/x];'j[|Ar-iVe,| >tViVh^] s;e^'"^(^-*'). (3.29) 

As an outcome of this article, we will be more precise in Section 8.2 about large deviations of filling 
fractions when the potential satisfies the stronger Hypotheses 1.1-1.4. 

3.4 Large deviation of Ln: distance (Lemma 3.5) 
3.4.1 Regularization of Ln 

We start by following an idea introduced by Mai'da and Maurel-Segala [MMS12, Proposition 3.2]. Let 
cntVn — > be two sequences of positive numbers. To any configuration of points Ai ^ . . . ^ Xm in 
A, we associate another configuration Ai, . . . , Aat by the formula: 

Ai = Ai, Aj+i = A,; + max(Ai+i - Ai, cr at), (3.30) 

It has the properties: 

ViT^j, |A, - Ajl > (TAT, |A, - Ajl «: |Ai-Aj|, |Ai - A,| (i - l)(TAr. (3.31) 

Let us denote Ljy = N^^ Sill '^a the new counting measure. Then, we define be the convolution 
of Ln with the uniform measure on [0, rj^aN]- 

We are going to compare the logarithmic energy of Lat to that of L^^, which has the advantage of 
having no atom. We may write by (3.31): 



^a{Ln) = // dLN{x)dLN(,y) In \x - y\ ^ // dLAr(a;)dLAr(y) In \x - y\ = SA(ijv) (3.32) 

and, if we denote lA^W are two independent random variables uniformly distributed on [0, 1], we find: 

SA(£jv)-SA(£iV) = / 'iLN{x)ALN{y)¥\\n{l + i^M<JN^^—^) 

Jx^y L V {x-y)J\ 

s; / dZAr(a::)dZAr(y) "^^^^ <c ry^^ 

Jx^y \x-y\ 
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thanks to the minimal distance between A^'s enforced in (3.30). Eventually, we compute with = 
//ln|a; - y\dn{x)dfi{y), 

I]a(Z5V)-S(Z^) = - // dL^^{x)dL^;,{y) ln\x - y\ 

JJ x—y 

= -^E[ln|,7jvaAr(Z^-Z^')l] = ^(^-31n(77^,a^^)). (3.33) 
Besides, if 6 > and (p : A ^ C is a 6-Lipschitz function with constant Kb{f), we have by (3.31): 

f d[L^ ZiVK-) ^(-)| < - l)'[-A^(l + ^m)]' < (Na^)' (3.34) 

A j^]^ V ) 

3.4.2 Large deviations of L'^ 

We would like to estimate the probability of large deviations of from the equilibrium measure 
Moq = fJ'Yq'^- We need first a lower bound on Z^'^ similar to that of [AG97] obtained by localizing the 
ordered eigenvalues at a distance of the quantiles of the equilibrium measure /it^''^, which are 
defined as: 

Af = inf{xeA, ^^^^^([-(X), x]) > (3.35) 



Since y is C^, d/ij^'''^ is continuous on the interior of its support, and diverge only at hard edges, where 
it blows at most like the inverse of a squareroot. Therefore, there exists a constant C > such that, 
for N large enough: 

\^f->^ti\^§2- (3-36) 
Then, since V is a fortiori on A compact, 

r ^ 



1 = 1 



N 

JV/3 Y^iV 



^ NlN-^^e-^^^ Yl \Xf-XffYle-^^"-^''^^'\ (3.37) 

lsii<j!^N i=l 



for some constant Ci > 0. Therefore, since: 



[[ ln|x-y|dMeT(x)dMeT(2/) < [ [ '"\n\x-y\d^^ll\x)d^^^^{y) 

< ^2 E ln|Af-Af| 

1 1 / \ 

< JpI.^^\^f-^f\+N^^{^)' (3-38) 

we find: 

ZnI^p > exp { ^ ( - C2 TV In iV - eI^^^]) } . (3.39) 

for some positive constant C'2 and with the energy introduced in (1.10). 
Now, let us denote 5jv(t) the event {D[Z^, ^J'^^'^j ^t}. We have: 
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and using the comparisons (3.32)-(3.34), we find, with the notations of Theorem 1.2: 



62 



N,/3 •'SN{t) 



N 



N,I3 



5jv(t) 



(s[LS,,M„j+/dL:iVW!^'''''(^)) Y\AXi, (3.41) 



where: 



3iV 



Rn = N^(TN Hl{V) + N^VN + ^ - 3iVhl(?7ArCrAr), 



(3.42) 



and the effective potential C/^'"^ was defined in (3.1). Since [/^'"^ is at least i-Lipschitz on A (and 
even 1-Lipschitz if all edges are soft), we find: 



7V;A 



r7 \ . . 
^ AJ a J5jv(t) 



N 



(3.43) 

We now use the lower bound (3.39) for the partition function, and the definition of the event 5Ar(t), 
in order to obtain: 



with: 



Rn = Rn + Ki/2(C/) N''l^a]i'' + ^ ln£(A). 



(3.44) 



(3.45) 



Indeed, since U^^'^ is nonnegative on A, we observed that the integral in bracket is bounded by the 
total length £(A) of the range of integration, which is here finite. We now choose: 



(3.46) 



which guarantee that i?jv e 0{Nh\N). Thus, there exists a positive constant C3 such that, for N 
large enough: 

M];;^[5^W]*:e^(^3^'"^-~^*^), (3.47) 

which concludes the proof of Proposition 3.5. We may rephrase this result by saying that the proba- 
bility of SN{t) becomes small for t larger than wn = -s/lC^ \nN/N. 



□ 



3.5 Large deviations for test functions 
3.5.1 Proof of Corollary 3.6 

Since (f is 6-Lipschitz, we can use the comparison (3.34) with = chosen in (3.46): 



< 



(6+ l)iV26 



(3.48) 



Then, we compute: 



d[Z^ - Meq](a;) f{x) 



ds(£^-Mcq)(s)^(-s) 
ds 



l¥'ll/2( / t4 K^'JV -Aioq)(s)|^) 



1/2 



(3.49) 
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and we recognize in the last factor the definition (3.22) of the pseudo-distance: 

/ d[Z^-^eq](a;)</?(a;)| < V2|</7|i/2 2)[Z^,//eq]. 

J ^ ' 
Corollary 3.6 then follows from this inequality combined with Lemma 3.5. 

3.5.2 Bounds on correlators and filling fractions (Proof of Corollary 3.7 and 3.8) 

If /i is a probability measure, let denote its Stieltjes transform. We have: 

[Wl„ - W^J{x) = [ A[Lm - McqKO MO, MO = (C) + = ^ 

Since tpx is 1-Lipschitz with constant ki(V'2:) = (i~^(a:, A), we have for N large enough: 



(3.50) 



(3.51) 



|[Wl„-W?„](x)U 



1 



We now focus on estimating W^u — W^^^^ . Since the support of is included in 

Ai/jv2 = {x e M, d{x,k) s^l/N^}, 



(3.52) 



(3.53) 



we have the freedom to replace ^/j* by any function which coincides with ?/;* on Aj^/jy^ • We also 
find: 

|[WZ"„ - WMeJ(^)| < V2 (|0f |i/2 + \4>i\ii2) (3.54) 
Let (j) : M ^ M be C""^ function which decays as 4){x) e 0(l/a;^) when oo. We observe: 

1 



I '/'1 1/2 



2 



|s||</.(s)|2ds = y^^|0'(s)|2rfs 
2 / In 16 - 61 </''(6)'/''(6)d6d6 = -2 /" In 



J^i^ </''(6)'/''(6)d6d6 



In 



M 



(3.55) 



where M > can be chosen arbitrarily. Let ax,h 6 f\h the point such that d{x,Ah) = \x — ax,h\- We 
claim that, for d{x,A) small enough, we can always choose (p^ and 0^ such that: 

9 



h=0 



1 



(3.56) 



This family (indexed by x) of functions is uniformly bounded by {g + 1)/^^ at oo, which is integrable 
at CO. Therefore, we can choose M in (3.55) independent of x so that: 



\€\l/2 *s 1 - ^2 ^^"^ (^Iw^) K'^*)'^^!)! I(<^*)'fe)|d6d6 



(3.57) 



If we plug the bound (3.56) in the right-hand side, the integral can be explicitly computed and we 
find a finite constant D > which depends only on A, such that: 



2 mnd{x^ 
I'^-li/^ ^ d2(x,A) 

when X approaches A. Combining (3.52)-(3.54)-(3.58) with Lemma 3.5, we find: 



(3.58) 



1 

N' 



w,ix) - wl: '\x) = A^];;^[WL„(a:) - w^„,(x)] 



Nd^{x,A) 



+ 2D 



\nN V'|lnd(a;, A)| 
d(x. A) 



(3.59) 
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If we restrict ourselves to x e C\A such that: 



for some constant D' , then: 



IniV ^/\liid{x,A)\ 
N d{x,A) 



^ W^{x)-Wt'H^) <(2^ + ^')\/^^^^P^- (3.61) 



iV 

Now let us consider the higher correlators. For any n ^ 2, WX''^ is the expectation value of some 
homogeneous polynomial of degree 1 in the quantities (Wln ^ W^/icq)(2^j) ^'^^ m]v'^ — yV^^^)(a;i)] . 

Accordingly, Lemma 3.5 yields: 



for some constant 7„ > 0, which depends only on A. This concludes the proof of Corollary 3.7. 
Similarly, to have a hand on filling fraction, we write: 



Nh - Ne.^h = N Jd[LN - /ieq](C) 1a, (0- (3.63) 

After replacing the function 1a, by a smooth function which assumes the value 1 on A/i, vanishes on 
Ah' for h' ^ h, and has compact support, we can apply Corollary 3.6 to deduce Corollary 3.8. □ 

4 Schwinger-Dyson equations for /3 ensembles 

Let A = Ufj^g Afi be a finite union of pairwise disjoint bounded segments, and 1^ be a function of 
A. Schwinger-Dyson equation for the initial model fJ,^'^ can be derived by integration by parts. Since 
the result is well-known (and has been reproved in [BGll]), we shall give them without proof. They 
can be written in several equivalent forms, and here we recast them in a way which is convenient for 
our analysis. We actually assume that V extends to a holomorphic function in a neighborhood of A, 
so that they can written in terms of contour integrals of correlators, and an extension to V harmonic 
will be mentioned in § 6.2. We introduce (arbitrarily for the moment) a partition dA = ((5A)_,_0(c5A)„ 
of the set of edges of the range of integration, and 

L{x) = II (a; -a), Li{x,Q = — , ^2(2^; 6) = ^-37 • (4-1) 

Theorem 4.1 Schwinger-Dyson equation at level 1. For any x 6 C\A, we have: 

2 



W2{x,x) + {Wi{x)y + (1- -)5,VFi(x) (4.2) 



And similarly, for higher correlators: 
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Theorem 4.2 Schwinger-Dyson equation at level n > 2. For any x,X2, . . . , x„ G C\A, if we denote 
I = [2, nj, we have: 



W„+i{x,x,xi) + 2 W^j^+i{x,xj)W„^\j\{x,xi\j) + (l- ^^d^Wn{x,xi) (4.3) 



A2 



The last line in (4.2) or (4.3) is a rational fraction in x, with poles at a e ((5A) + , whose coefhcients 
are linear combination of moments of A^. 

As a matter of fact, if AT e [0, so that \N\ = X;^=i < the correlators in the model with 
fixed filling fractions M/v e^jv/Af /3 satisfy the same equations. Indeed, in the process of integration by 
parts, one does not make use of the information about the location of the A^'s. By linearity, in a model 
TV where is random, the partition function Z = ]E[z]^'^y^^] and the correlators Wn{xi, . . . ,a;„) = 
^[^N/N.n{^i: ■ ■ ■ j^n)] Satisfy the same equations. The initial model /i]^'^ and the model with fixed 
filling fractions are just special cases of the model with random filling fractions. 

When g > 1, we denote ^ = {Ah)i^h!ig a family of contours surrounding in C\A, and introduce 
the vector-valued linear operator: 



(4.4) 



on the space of holomorphic functions in C\A. For any me [1, We denote: 

) = /:^'"[l^„(xi,...,x„_™,.)], (4.5) 

which means that we integrate the remaining m variables on ^-cycles. By definition, if we denote 
(e/i)i^/i<;g the canonical basis of C^, 

n n-m 1 1 ™ 

W„|™(xi,...,a;„_™) = ^ Mp,t,J n ^/^^ n Tr — -r (X)e;,,. (4.6) 

ls;/ll,...,h,„^g i = Tl-T?l+l j = l J J = l 

In particular, W„|„ is the tensor of n-th order cumulants of the numbers of A's in the segment 
Ah- We take as convention Wn\o{xi, ■ ■ ■ , Xn) = Wn{xi, . . . , a;„). Here, fJ^^'^ ^ denotes the probability 
measure in a model with N A's and random filling fractions, li e e £g and if we specialize to the model 
with fixed filling fraction e, we have: 

Wn\mixi, Xn-ni) = 5n,lSm,l Ne. (4.7) 

In general, we deduce from Theorem 4.2: 

Corollary 4.3 For any n ^ I and me |0,7i — 1], for any x, X2t ■ ■ ■ , Xn-m G C\A, if we denote 
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/ = [2, n — to], we have: 

Wn + i\,ri{x,X,Xl) + Yj ["\]'^\J\ + l + ni'\ni'{x,Xj)®Wn-\.]\-m'[m-m'{x,Xj\j) (4.8) 



ael^A), ^ - « "'"'^"-^'"^ + /? I-Ja 2i7r L(x) (X - a(x, - 
2x r i2(x;^,C) , ir dad6 i2(x;6,6) 



^V-Vt"^ -^W~ " Fa^ W L(x) X,) 

+ Yj ( ™, )^./|+i+"i'|ni'(Ci7a;,/) ® W„_|.7|_„.|„_„'(C2,a;/\j)^ 



Proof. Straightforward from (4.3), once we notice that the integrals over a closed cycle of a total 
derivative or a holomorphic integrand in neighborhoods of A in C give a zero contribution. □ 

5 Fixed filling fractions: expansion of correlators 

5.1 Norms on analytic functions and assumptions 

In this Section, we analyze the Schwinger-Dyson equations in the model with random filling fractions 
(which contains the model with fixed filling fractions as special case) and the following assumptions: 

Hypothesis 5.1 

• A is a disjoint finite union of bounded segments Ah = [a^ , a^] . 

• (Real-analyticity) V : A — > M extends as a holomorphic function in a neighborhood U c C of A. 

• (1/N expansion for the potential) There exists a sequence (V^''^)k^o of holomorphic functions 
in U and constants (w''^^)fe>0) ■so that, for any K ^ 0; 



sup 



K 



k=0 



^(0 - 2 ^"''^^''^(C) s; w*^^ A^"(^+^^ (5.1) 



» (g + 1-cut regime) ^\x) = liuiN^rc Wi{x) exists, is uniform for x in any compact 
of C\A, and extends to a holomorphic function on C\S, where S is a disjoint finite union of 
segments Sh = [Q!^,a^] c: A/j. 

• (Offcriticality) y{x) = '-^ ^ — Wi '^\x) takes the form: 



y{x) = S{x) Yl^{x-a+)Pi{x-a-r'., (5.2) 



h=0 



where S does not vanish on A, are all pairwise distinct, and pj^ = 1 z/a^ e dA, and = — 1 
else. 
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(1/N expansion for correlations of filling fractions) For any n ^ 1. there exist a sequence 
(W^^2)'c>n-2 of elements of (C^)®", and positive constants {w^l'^j^)k^n-2, so that, for any K > 
— 1, we have: 

k=n-2 



W„|„- y 7V-^VF<f* 7V-(^+i) u;*;"* (5.3) 

n\n / ] n\n n\n ^ 



We say that Hypothesis 5.1 holds up to o{N ^) if we only have a expansion of the potential at 
least up to o{N^^), and an asymptotic expansion for correlations of filling fractions up to o(A^^^'^^^^). 

Remark 5.2 In the model with fixed filling fractions, this last point is automatically satisfied since 
Wn\n is given by (4.7). Then, ^\x) is the Stieltjes transform of the equilibrium measure deter- 
mined by Theorem 1.2, and Hypothesis 5.1 then constrains the choice of V and e. 

We fix once for all a neighborhood U of A so that S'^^(O) n U = 0, and contours A. = {Ah)i<ih!^g 
surrounding A/i in U. 

Definition 5.3 If S > 0, we introduce the norm \\ ■ \\g on the space 'Hml^...,mn{^) of holomorphic 
functions on (C\A)" which behave like 0(l/a;™') when Xi —* co: 

\\f\\s= sup |/(a;i, . . . ,a;„)| = sup |/(a;i, . . . , x„)|, (5.4) 

d{x,A)^S d{xi,A) = 5 

the last equality following from the maximum principle. If 2, we denote Hm^ = 'Hm}...,m. 

From Cauchy residue formula, we have a naive bound on the derivatives of a function / e in 
terms of / itself: 

11^-/(^)11. <^^ll/ll./2- (5-5) 
Definition 5.4 We fix once for all a sequence 5m of positive numbers, so that: 

lim^(^)^^^^0 (5.6) 

If f & Hm'*(A) is a sequence of functions indexed by N, we will denote f 6 0{Rn{6)) when, for any 
e > 0, there exists a constant C(e) > independent of d and N, so that: 

WfWs ^ C{s)6-' Rn{S) (5.7) 
for N large enough and 6 small enough but larger than Sn ■ 

This choice will be justified at the end of § 5.3.1, and we can simplify the condition by taking Sjy of 
order A^^^/^+s for some e > arbitrarily small. We notice it also guarantees the assumptions d{x,A) 
not much smaller than (N \n N)~^^'^ as it appears in Corollary 3.7. 

Definition 5.5 If X is an element o/(C^)®" , we define its norm as: 

\X\= |^/.„...,/.J. (5.8) 

l!ihi,...,h„!ig 
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To perform the asymptotic analysis to all order, we need a rough a priori estimate on the correlators. 
We have established (actually under weaker assumptions on the potential) in § 3.3 that: 

(Wi - ^ T^i*"^^) e 0{^/N\nN5~^) (5.9) 
and for any n ^ 2 and m 6 [0, n]: 

W„|„ e 0((iVlnAf)"/2^-("-™)«) (5.10) 

with exponent 6=1. 

Our goal in this section is to establish under those assumptions Proposition 5.5 below about the 
1/iV expansion for the correlators. We are going to recast the Schwinger-Dyson equations in a form 
which makes the asymptotic analysis easier. We already notice that it is convenient to choose 

{dA)+ = {ale{dA), pl = ±l}, (5.11) 

as bipartition of dA to write down the Schwinger-Dyson equation, since the terms involving da In Z 
and daWn-i for a e ((5A)+ will be exponentially small according to Corollary 3.3. If a = a* , we denote 
a{a) = al- 

5.2 Some relevant linear operators 
5.2.1 The operator K. 

We introduce the following linear operator defined on the space T-6^\A): 

Km = 2wt\.)m - i m^^mqm ^ p{-ii(,;^)i;(e), (5.12) 

Lyx) Jf^ 2i7r L a; — ^ J 

where: 

p{-i} (a;; 6 = ^ ^ ^L2{x; e, ry) W^'^ (77) (5.13) 

We remind that L(x) = Y\ae(dA) (-^ ^ a{a)) and L2 was defined in (4,1). Notice that ^\x) - 1/x 
when a; ^ 00, and P^^^^ {x,S,) is a polynomial in two variables, of maximal total degree |((?A)_| — 2. 
Hence: 

IC : H'i\A) ~^h['\a). (5.14) 

Notice also that: 

,„.fiq)M-,^!-«„.5wy||. (5,15, 

where L(x) = na€(c5A)+(^ ^ &n.d by offcriticality assumption the zeroes of S are away from A. 



Let us define (t{x) = ^^jYlae{dA)(^ ^ a(a)) = y L(x)L(x), so that ^||jy = §|fy- We may rewrite: 

JCf{x) = -2y{x)f{x) + ^/^, (5.16) 

L{x) 

where: 

J A 2i7r L X — t J 
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For any / e 'H'^\a), Qf is analytic in C\A, with singularities only where has singularities, 

in particular it is holomorphic in the neighborhood of A. It is clear that Im/C c 'H^^'*(A). Let (p elmJC, 
and / e 'H''2\^) such that </? = /C/. We can write: 

a{x)f{x) = Res ^(0 /(C) = VX^:) - (L ^ "f^^^K (5.18) 

where: 

V'(x) = - Res -^a(C)/(a. (5.19) 

Since /(x) e 0(l/x^), ip{x) is a polynomial in x of degree at most <? — 1. We then pursue the 
computation: 



jTi 2i7r C - X 2S'(C) 

using the fact that 5 has no zeroes on A. Let us denote G ■ Im/C 'H'^\a) the linear operator 
defined by: 

mix) = J-(f^^I^ ^(0. (5.21) 
a[x) Jf, 2m £,-x 2S[Q 

One also obtains: 

/(^) = ^ + (^o/C)[/](x). (5.22) 



It was first observed in [Akc96] that ^^y)^)^ defines a holomorphic 1-form on the Riemann surface 



5.2.2 The extended operator JC and its inverse 

j{x) 

E : (7^ = riae(5A)('^ ^ a{a)). The space -ff^(S) of holomorphic 1-forms on E has dimension g if all 
a{a) are pairwise distinct (which is the case by offcriticality) and the number of cuts is {g + 1) and. 
So, if g > 1, /C is not invertible. But we can define an extended operator: 

£ : -H^^^A) — > Im/C x C 

/ ^ {Kf,C^[fdx]). (5.23) 

Since { "^-^T^r] are independent, they form a basis of H^{Y?). On the other hand, the family 

of linear forms: 

CA=i$ ,---,S ) (5.24) 



are independent, hence they determine a unique basis Wh{x) = ^''J^^^^ e H^{T.) so that: 

m,h'ell,gl (j) Wh'{x) =Sh.h'- (5.25) 

Cj\, thus induces a linear isomorphism of the space of H^(Yi). Its inverse can be written: 

a 

C^iw] = ^ WhWh{x) (5.26) 

h=l 
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We deduce that JC is an isomorphism, its inverse being given by: 

IC '-[ip,w]{x) = +Gip{x), (5.27) 

ax 

where G is defined in (5.21). We will use the notation JC^ip = K.^^'up^w]. The continuity of this 
inverse operator is the key ingredient of our method: 

Lemma 5.1 Im/C is closed in TL^^A), and there exists s constant fcr(A), fcp(A) > such that: 

V((^,w) 6 Im/C X ||£^Vllr<^r(A)||(y5||r + A:i.(A)|i(;| (5.28) 

□ 

Remark 1. If one is interested in controlling the large N expansion of the correlators explicitly in 
terms of the distance of Xi's to A, it is useful to give an explicit bound on the norm of IC^^. For this 
purpose, let Sq > he small enough but fixed once for all, and we move the contour in (5.21) to a 
contour close staying at distance larger than 5q from A. If we choose now a point x so that d{x, A) < 77, 
we can write: 

C ( )= -^cf ^ ^ ^ + ^1 K ^ ^ (K oq) 

^ 2S{x)a{x) <j{x) LaA)=So 2i7r 2S{0 x C <j{x) LiiA)=So 2i7r 23(0 x-^ ^ ' ^ 

Hence, there exist constants C, C" > depending only on the position of the pairwise disjoint segments 
Ah such that, for any S > smaller than So/2: 

\\g^\\, ^ {CD,{5) + C) 5-^1^ + r ii^ii^^ (5.30) 

We set: 

D,{S) = sup ^ (5.31) 



For S small enough but fixed, Dc{5) blows up when the parameters of the model are tuned to achieve 



-1 



a critical point, i.e. it measures the distance to criticality. Besides, we have for the operator 
written in (5.2G): 

||/:,nHII.^=^^^^*|H, (5.32) 

inid(4,A)=5|CT(a;j| 

and the denominator behaves like J"^/^ when 5 ~* Q. We then deduce from (5.27) the existence of a 
constant C" > so that: 

||£„Vll5 < {CD,{5) + C')5-'^M\s + 5-'^ \wl (5.33) 

with exponent k = 1/2. 

Remark 2. From the expression (5.27) for the inverse, we observe that, if Lp is holomorphic in C\S, 
so is ^lo^yj for any w e C^, in other words £~-^(Im/C n 'H^i\S)) c 'h':^\S). 

5.2.3 Other linear operators 

Some other linear operators appear naturally in the Schwinger-Dyson equation. We collect them 
below. Let us first define: 

A_iT4^i(x) = N-^Wi{x) -W[~^\x), (5.34) 

A_iP(a;;0 = (f) ^ 2L2{x; ^,v)A_,Wiir,), (5.35) 
/a 2i7r 

AoV{x) = V{x)-V^°\x). (5.36) 
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Let also be two holomorphic functions in U. We define: 



'"i 








C2 




-^h[^\a) 


C2f{x) = ^ 






^Hf\A) 


Mx'fix) = { 






-^n[^\A) 


■^hiMfix) = 


A/C 


h[^\a) 


-^n[^\A) 


(A/C)/(x) = 



(2i^)2 L{x) 
1 / de 



L(a;) 2i7r V a; — ^ 



(5.37) 



All those operators are continuous for appropriate norms, since we have the bounds, for small 
enough but fixed, and 8 small enough: 



ll(A/C)/| 



^^llfll 



^ ll/^rll? 11/11. + CMML±I«||;||^,^, 



(||(Aoy)'ir^ + 2||A_iiyi||,)||/||,+ 
ll £(AoV^yHg' + HA_iP||'^. 



2 



2C 



11/11, 



ll-' 115/2 



+ C 



for any / in the domain of definition of the corresponding operator, and: 

C = £(A)/7r+(g+l), I?l(<5)= inf \L{x) 



(5.38) 



(5.39) 



If all edges are soft, Dl{S) = 1, whereas if there exist at least one hard edge, Dl{5) scales like 6 when 
,5^0. 



5.3 Recursive expansion of the correlators 
5.3.1 Rewriting Schwinger-Dyson equations 

For n 5= 2 and me |0,n — 1], we can organize the Schwinger-Dyson equation of Corollary 4.2 as 
follows: 



/C + A/C + l(l-|,.. 



(5.40) 
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where: 

Ai+i\,n{x;xi) = N^^{C2 - id)Wn+l\„i{x,X,Xi), 

{x,Xj) (^Wn-\J\ (a;,x/\j)|, 

(J,m')?^(0.O),(/,m) 

2 

C„_i|m(a;;x/) = --^^ A^2;,W„_i|„(a;,x/\{i|), 

-D„_iim(a;;x/) = 1^ Yi ^i^ri-i\rnixi). (5.41) 

This equation can be rewritten, with the notation IC^Lp = IC^^ [ip, w] and by definition (4.5) of W„j„, 

''^n\rn{Xy Xj) = A]!^^^^^^^^^^ |^^„^x|,„ (x, X/) + B„|„j (x, ) + Cn—l\m{x,Xj) + Z)„_]^|^(x, X/) 

- {AIC)[W,,\^{x,xi)] - ^(l - |)((?. + Ci)Wr,\m{x,Xl)], (5.42) 
where it is understood that the operators ali act on the first variable (or the two first variables for 

For n = 1 and to = 0, we have almost the same equation; with the notation of (5.34), and in view 
of (4.2), 

A,,W,ix) = ^-A\.w4^^^^-j^{l-^yS. + C,)W,ix) 

-^(AoVYfiWiix) - (A_iW^i(x))'l, (5.43) 



where: 



A_iVFi|i = Ca[A-iWi] = N-^ £AWi] ~ Ca\w[ (5.44) 



is the first correction to the expected filling fractions. 

We would like A/C to be negligible compared to K in (5.40), that is 

||£^(A/C/)||, « 11/11, (5.45) 

This can be controlled thanks to (5.33) and (5.38). From our assumptions, we know that: 

||A_iP||;^. e 0(1), ||AoV^||;^ eO(l/7V). (5.46) 

Taking into account those estimates in (5.38), we observe that it will be possible independently of the 

— 1/2 

nature of the edges provided b is restricted to be larger than Jat such that limTv^ac ^ ^ 



1-2//3 

independently of /3 if: 



0, and limjv-joo = 0- Given the a priori bound in Corollary 3.7 on A_iW^i, this can be realized 

"if: 

/in^Vre^^^^^^ (5-47) 

Since we consider here a fixed, off-critical potential, Dc{S) remains bounded, and since k = 1/2, this 
condition is equivalent to: 

/lniV\i/3 ln^/3,5^ 
lim , = 5.48 

N^oo \ N J 6n 

It justifies the introduction of the sequence J^r in Definition 5.4. Then, by similar arguments, ||D„|„j||, 
will always be exponentially small when N is large provided S ^ Sn, and thus negligible in front of 
the other terms. 
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5.3.2 Initialization and order of magnitude of Wn 

We remind 9=1 and k = 1/2 here. The goal of this section is to prove the following bounds. 
Proposition 5.2 There exists a function Wi'^^ e 'H'^\S) independent of N so that: 



Wi = NW^ ^* + + AoWi, A, 



w!°Hx) = K-^,,, {[- (l - |)(^- + A) -^^ivi-yy,o\wt'^^ (x) (5.50) 



It is given by: 

P 

Proposition 5.3 For any 2, there exists a function Wn^ e 'h!:^\S) so that: 

t^„l™ = ^'^"(W^ir„ + A„_2W„i^) (5.51) 

where: 

A„_2W^„ 6 0(iV-i/2(ln7V)2"-3/2 (^^^)'""' ^-e(„-™)-(«+e)(3„-3)^ (5 52) 

Prior to those results, we are going to prove the following bound: 

Lemma 5.4 Denote r* = 3n — 4. For any integers n 5= 2 and r > smc/i i/iat r ^ r*, we have: 

I f iV^ (In TV) ^ f r '^("-™)-(^+^)'^) (5.53) 

V \DjXd)J ) 



Proof. The a priori control of correlators (3.28) provides the result for r = 0. Let s be an integer, 
and assume the result is true for any r e [0, s]. Let n be such that s + 1 r* = 3n — 4. We consider 
(5.42) which gives W„|j„ in terms of W„+i[m f^nd Wn'\m' for < n or n' = n but m' > m, and we 
exploit the control (5.33) on the inverse of /C. We obtain the following bound on the A-term: 



^ 0(A^^(lnA.)^ + ||T^„j„^^,|l^<5-V.), 

(5.54) 

and we now argue that it will be the worse estimate among all other terms, given that S ^ S^. The 
_B-term involves linear combinations of Wj+i|m' ® ^n-j|m-m'- Notice that: 

s — r*^-^ ^ rj^j — r*^i = r*_j (5.55) 

Therefore, we can use the recursion hypothesis with r = r*_^^,_^^ = 3(j + m') — 1 to bound Wj^jn'+i\m' ^ 
and with r = s — r*j^^^,^^ to bound Wn_j_„i'\m-m' ^ and we find: 



-l^^JB^\m) e 0(iV^ (IniV)^ (^)"' + \\W^\.^,,\l 5^) 

(5.56) 

which is of the same order as (5.54). The C-term involves W„_i|„j. If s ^ '^^ •^^'r ^^e the 

recursion hypothesis at r = s to bound it as 
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The prefactor makes this term neghgible comparing to (5.54). If s > r*_^, we can only use the 
recursion hypothesis for r = r*_^, and find the bound: 

+ l|W'„|,„+i||,r«), (5.58) 



and thus still negligible compared to the A-term when N ^ co and S ^ 0. We can conclude by 
recursion from to = n to m = 0, taking into account the assumed expansion (5.3) for W,j|„. □ 

Proof of Proposition 5.2. Lemma 5.4 for r = 1 gives the bound: 

W2eO(VNh^^^6-^''+-A (5.59) 

Then, in (5.43), we find that A2 e O ( jv_d1V) ^ ' ^^'^ already know that |A_iVFi| goes to zero 
for uniform convergence in any compact subset of C\A, so that ^ (A_iT4^i)^^5 is neglectable 

with respect to A_iVFi at list for (5 not going to zero with N . By continuity of we deduce 

that A''A_iWi(a;) has a limit wf'^ for convergence in any compact subset of C\A, given by (5.50). 
Reminding Remark 2 page 28, this limit belongs to 'H2^'*(S), and given the behavior of w\ at the 
edges, we have a naive bound: 

(c5, + /:i)W^reO(£^), (5.60) 

and we have already argued at the end of § 5.3.1 that ^N~'^(dx + >Ci)(A_i Wi) ||^ was negligible 
compared to ||A_iWi 11^/2 provided & 5= Moreover, Therefore, the worse estimate on the error 
AqW^i is given by the term A2. Taking into account the effect of Kr^ given in (5.33), we find: 

which is the desired result. □ 

Proof of Proposition 5.3 We already know the result for n = 1. Let n 5= 2, m e [0,n — 1], and 
assume the result holds for all v! e — 1] and to' e [TO,ri]. We want to use (5.40) once more to 
compute Wii[,„. Applying Lemma 5.4 to r = 3n — 4 for n 5= 2, we find: 



= + A«-2W„u„) with: 



whereas the recursion hypothesis implies that B„|„j and Cn-\\m are of order 0(N^ "). Hence, 



,7Q/ 



The error term lS.rL-2^n\m receives contribution either from errors A„'_2W„/|m/ appearing in S„[^ 
and C„_i|to, which we already know how to bound. And the restriction that limjv^oo A^<5^ = +Q0 
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guarantees that they are negUgible in front of An+i\m computed in (5.62). Hence, we deduce by 
substracting W^^^^ and inverting /C that: 

A.-.T^H™^o(^/(1^^(^)^^^^ (5.64) 
which is the desired result at step (n,m). We conclude by recursion. □ 
5.4 Recursive expansion of the correlators 

Proposition 5.5 For any ko > 0, we have for any n > 1; 

W„l„(xi,...,x„) = J] N-''wj^l{xi,...,Xn)+N''^"{Ak,W^^^){x,,...,x^), (5.65) 

k=n-2 

where: 



(i) for any n > 1 and any k e [O,...,fcol; ^'^^ limit when iV — > co in 'H'^ '"■'(S) for 

pointwise convergence in any compact of (C\A)"^™, and: 

VF<^^ e of(ln7V)^ f#^)"^''<5-("-'")«-("+2'=)(«+«)y (5.66) 
"I"* v \D]^yo)J J 

{ii) for any 1, AkoWn\m e "^\a) and: 

Proof. The case fcg = follows from § 5.3.2, and we prove the general case by recursion on fco, which 
can be seen as the continuation of the proof of Proposition 5.3. Assume the result holds for some 
ko > 0. Let us decompose: 

ko + 2 

V = N-'' V^''^ + N-^''°+^^ Ak,+2V. (5.68) 

fc=0 

We already know that the loop equations are satisfied up to order N^~'^°. We can decompose the 
remainder as: 

^-iko-i) |/c + A/C + ^(l - |)5,}Afc,W„|„(x,x,) = N-'''>{Ei^^°Jix;xi) + H^^'(x;x,)). (5.69) 

It is understood that all linear operators appearing here (and defined in § 5.2) act on the variable(s) 
X. We have set: 

= (^■2-^<i)[wi%:^\x,x,xj)] 
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and 

>{fco} 



J\—m'\7n—m' 



;=-i P iel 

2 ■r-i L{a) 
-"5 >, (5aW;i-iim(a;/ . 

aG(f?A) + 



Let us denote: 



Thanks to the recursion hypothesis and the bound on the norm of the inverse of /C from (5.33), we 
find: 

iCo'Rn\,neO{wi'+'^iN,S)). (5.72) 

This bound arise from the three first hnes. Indeed, the two last terms are negUgible compared to 
wl^^^^ as noticed in the proof of Lemma 5.4, and the term involving Ai+iV is of order (In A^)^/A^ for 
some p, with a dependence in S which is less divergent than that appearing of Wn^^^K Therefore, we 
deduce by recursion on m from m = n to m = that NAf^gWn\m converges to: 

pointwise and uniformly on any compact of (C\A)"^™. Besides, the estimate (5.72) yields the bound 
(5.67) for the error, while the recursion hypothesis combined with (5.73) leads to (5.6G). □ 
This proves the first part of Theorem 1.3 for real-analytic potentials (i.e. the stronger Hypoth- 
esis 1.3 instead of 1.4). For given n and k, the bound on the error A^Wn depend only on a finite 
number of constants v^'^ \ w^, ' appearing in Hypotheses 5.1. 

5.5 Central limit theorem 

With Proposition 5.2 at our disposal, we can already establish a central limit theorem for linear 
statistics of analytic functions in the fixed filling fraction model. It will be refined for non-analytic 
but smooth enough functions in § 6.1. 

Proposition 5.6 Assume the result of Proposition 5.2. Let : A — > M extending to a holomorphic 
Junction in a neighborhood of S . Then: 

N 



exp ( ^(A,))] = exp (iVLM + M[ip\ + ^ Q[^, ^] + o(l)) , (5.74) 

i— 1 



where: 

L[^] = ^ m wt'HO, M[h] = ^(0 Wl'\0 (5.75) 



J A 2i7r 7^ 2i7r 

Wi^^ has been introduced in (5.49), and Q is a quadratic form given in (5.78) or (5.79) below. 
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Proof. Let us define Vt = V — ip. Since the equilibrium measure is the same for Vt and V, we 
still have the result of Proposition 5.2 for the model with potential Vt for any t e [0, 1], with uniform 
errors. We can thus write: 



1 v-A 



= £di^||^(e)[^<'^<-'U<'^*°>(e)]+o(l) (5.76) 



As already pointed out, W^''^"'^ = W^'^°\ and from (5.49): 



^v.io} ^ ^vm 2t OA/-,, o)[<^<->] (5.77) 

P i|i 



Hence (5.74), with: 



Qb, (^] = - 1 ^ ^ ^{0 (£^,0, o K',o) [Wp^-'h iO (5.78) 



If we restrict to the model with fixed filling fraction, it can be simplified to: 

Qi^^f] = (5.79) 

where H^2^'^"' has been introduced in (5.51) and we recall W^^^ = for the model with fixed filling 
fractions. From the proof of Proposition 5.2, we observe that the o(l) in (5.74) is uniform in h such 
that sup^j^p^ l¥'(C)l is bounded by a fixed constant. □ 

In other words, the random variable $ = YiiLi fi^i) ~ ^[f] converges almost surely to a Gaussian 
variable with mean M[(/3] and variance (5[(/9, (^]. This is a generalization of the central limit theorem 
already known in the one-cut regime [Joh98, BGll]. A similar result was recently obtained in [Slicl2]. 
In the next Section, we are going to extend it to holomorphic h which could be complex-valued on A 
(Proposition G.3). In general, to establish the central limit theorem, one could be tempted to use the 
definition of the correlators: 



N „ n 

exp(2i^(A,))]| „ = n 



(5.80) 



then represent Gpfit) = m]v'^ [c**] by its Taylor expansion up to i = 1, and use the result of 
Proposition 5.2 that e 0{N^~'") to conclude. However, for any fixed N, GN{t) is analytic in the 
domain of the complex plane where At]^'^ ^ [g**] does not vanish, and it is not obvious that for N large 
enough (although it will turn out to be true) that this does not happen for some <o 6 C with |io| < 1, 
i.e. that the Taylor series converges in the appropriate domain. 



6 Fixed filling fraction: refined results 

In this section, we show how the asymptotic expansion of multilinear statistics for non-analytic test 
functions can be deduced from our results, thanks to their explicit dependence on the distance of the 
variables x (appearing in the correlators) to A. We also show how to extend our results to the case of 
harmonic potentials, and potentials containing a complex-valued term of order 0(1/A^). The latter is 
performed by using fine properties of analytic functions (the two-constants theorem) as was recently 
proposed in [Shcl2]. 
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6.1 Multilinear statistics for non-analytic test functions 



Our methods establish a control on n-point correlators Wn{xi, . . . , a;„), depending on how xi, . . . , a;„ 
approach the range of integration A. We now argue that it gives a control n-linear statistics for 
test functions with regularity lower than analytic. If s is a finite-dimensional vector, we denote 
|s|i = I],k.|. 

Lemma 6.1 Let fn be a holomorphic function defined in a neighborhood of /V^ in (C\A)". Assume 
there exists C,r, and rj 6 (0, 1) small enough, such that 



\fS^j], sup |/„(Ci,...,60K 



Then, there exists a constant C" so that, for any s satisfying \s\i e [r,r/r]\, we have: 



(6.1) 



(6.2) 



Proof. For 5 small enough but larger than rj, let C{S) be the contour surrounding A such that 
d(f , A) = S for any ^ e C{S). When A has (5+I) connected components, its length is 2(£(A) + (g+l)7r(5). 
For any s e M, we find: 



i n 



2i7r 



sS C 



C"(5) 

'£(A) 



^e-'^^«^/„(Ci,...,?„) 

ZITT 



(^^ + (,g+l)^)%N^^--'^^ 



(6.3) 



We now optimize this inequality keeping |s|i large in mind, by choosing 5 = r/|s|i, which leads to the 
desired result. □ 

Corollary 6.2 Let Lp : M" — > C &e a continuous function with compact support, so that its Fourier 
transform satisfies: 

^{s)^o{\s\--), \s\^^ (6.4) 
Then, for any integer ko such that r^l + n + 26+ ^(k + 0), we have an expansion of the form: 



n N ko 



(6.5) 



j = l = l 



k=n~2 



Proof. Let r] > 0, and define a function 1^9^ by its Fourier transform <fri{s) = e ■?(*)• It is analytic 
in the strip e C, |Ini^| < rj}, and we may write: 



V: 



P n p n 

^ I(n/''--""'«">)l(n 



2i7r 



(6.6) 
)w^n(ei,...,Cn) 



We may insert the large N expansion of the correlators established in Proposition 5.5: 

iy„(6,...,U= X] iV-'W^i'*(ei,---,Cn)+iV-'=«Afc„W^„(Cl,...,Cn) 

k=n-2 



(6.7) 
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where: 



We may pass to the limit 7y — > in (6.6) when the integrand in the right-hand side is integrable near 
|s|i = 00. It constrains the allowed behavior for (f{s) at \s\ co. The worse behavior at |s|i = oo 
comes from the error term AkgWu- Lemma 6.1 implies that, for any e > 0, there exists a constant 
Ce > such that: 



(6.9) 

Assume now that <p{s) e o(|s|^''). Then, integrability at |s|i in (6.6) requires: 

ne+{n + 2ka){K + 6)- n{r + e) < -n (6.10) 
In other words, performing an expansion up to o{N^^") if the regularity exponent r satisfies: 

rl^l + K + 29+—{K + 9). (6.11) 
n 

□ 

6.2 Extension to harmonic potentials 

The main use of the assumption that V is analytic came from the representation (1.6) of n-linear 
statistics described by a holomorphic function, in terms of contour integrals of the n-point correlator. 
If If is holomorphic in a neighborhood of A, its complex conjugate Tp is antiholomorphic, and we can 
also represent: 



N 

V;i 

1=1 



r A 

1 [ 2 Rao] =j)^^^{x)W,{x) (6.12) 



In this paragraph, we explain how to use a weaker set of assumptions than Hypothesis 1.3 , where 
" analyticity" and "1/iV expansion of the potential" are weakened as follows. 

Hypothesis 6.1 • (Harmonicity) V : A ^ R can be decomposed V = Vi + V2, where Vi,V2 
extends to holomorphic functions in a neighborhood U of A. 



(1/N expansion of the potential) For j = 1,2, there exists a sequence of holomorphic functions 
(vj'^^)fe^o cind constants (wj'^^)fc so that, for any K ^ Q: 

K 



sup 



- 1] N-"" Vf '(C) v]""^ iV-(^+i) (6.13) 



fc=0 



In other words, we only assume V to be harmonic. "Analyticity" corresponds to the special case 
V2 = 0. The main difference lies in the representation (6.12) of expectation values of antiholomorphic 
statistics, which come into play at various stages, but do not affect the reasoning. Below chronologi- 
cally Section 5, we enumerate below the small changes to take into account. 

In § 4, in the Schwinger-Dyson equations (Theorem 4.2 and 4.2), we encounter a term: 
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It is now equal to: 



' ^''mmwj(..,)-^S^mmwj(..,i (6,15) 



Remark that (6.14) or (6.15) still defines a holomorphic function of x in C\A. The second line in 
Corollary 4.3 has to be modified similarly. In § 5.2, we can define the operator /C by (5.16) with Q{x) 
now given by: 

^ r d^ m{vi'^y{o-L{x){vi'^y{x) ^^^^ 

/a 2i7r (-X 

^ideL(0(vf)-(0-L(.)(vf >)-(.) 
/a 2i7r 4 - X 

It is still a holomorphic function of x in a neighborhood of A, thus it disappears in the computation 
leading to formula 5.22 for the inverse of /C, which still holds. In § 5.2.3, the expression (5.37) for the 
operator A/C used in (5.40) should be replaced by: 



(A/C)/(x) = 2A_iW^i(x)/(x) + l(l-|)£i/(x) 



-A/'(AoVi)'.A_iP(^;.)[/](a;) - A/'(AoV2)',o[/](^), (6-17) 

and the bound (5.38) still holds, where vq is replaced by vi^ + V2,o introduced in (6.13). In § 5.3.1- 
5.4, all occurences of A/V',o[/](a;) should be replaced by ■^f(Viy,o[f]{^) + -^(V2)',o[/](^) (^-nd similarly 
for A/(AfcV)',o or J\fi^y{k}y q). The key remark is that the terms where V2 appear involve complex 
conjugates of contour integrals of the type g{0^ri''\^jXi) or g{^) AkWn{^,xi) where g is some 
holomorphic function in a neighborhood of A. Their norm can be controlled in terms of the norms 
of Wn'^^ or AkWn on contours F, as were the terms involving Vi, so the recursive control of errors in 
the 1/N expansion of correlators for the fixed filling fraction model is still valid, leading to the first 
part of Theorem 1.3, and to the central limit theorem (Proposition 5.6) for harmonic potentials in a 
neighborhood of A, which are still real-valued on A. 

6.3 Complex perturbations of the potential 

Proposition 6.3 The central limit theorem (5.74) holds for if : A — > C, which can be decomposed as 
(p = ip + where tpi, ip2 are holomorphic functions in a neighborhood of A. 

Proof. We present the proof for 93 = t /, where t e C and / : A ^ M extends to a holomorphic 
function in a neighborhood of A. Indeed, the case of / : A ^ M which can be decomposed as 
/ = /i + /2 with /i,/2 extending to holomorphic functions in a neighborhood of A, can be treated 
similarly with the modifications pointed out in § 6.2. Then, ii ip : A — > C can be decomposed 
'AS ip = (fii + p2 with (^1,(^2 holomorphic, we may decompose further pj = ip^^ + ipj, then write 
V = V~ ^(^f + ^f) and f={ipi- p^i), and: 

N N N 



ViA 



exp (2 MA,))] = A^P,t^[exp (2(v'f + V?))] 4;t^[exp i/(A.))]- (6-18) 



The first factor can be treated with the initial central limit theorem (Proposition 5.6), while an 
equivalent of the second factor for large N will be deduced from the following proof applied to the 
potential V. 
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This proof is inspired from that of [Shcl2, Lemma 1]. From Theorem 1.3 appHed to V up to 
0(1), we introduce wij"] for (n, fc) = (1, -1), (2, 0), (1, 0) (see (5.51)-(5.49)). If t e M, the central hmit 
theorem (Proposition 5.6) apphed to (p = t f imphes: 

N ,2 

I^N!''A{T,tf(^^))]=GNim + RNit)), G^(i)=exp(7VtL[/]+0/[/] + -Q[/,/]), (6.19) 
where sup^^|-_7^^ rp^^ \R]\[{t) \ ^ C{Tq) rjj^ and limjv^oo W — 0. Let Tq > 0, and introduce the function: 

RN{t) = \ RN{t). (6.20) 
For any fixed N , it is an entire function of t, and by construction 

sup |E7v(i)|s;l. (6.21) 

te[-To,To] 

Besides, for any t e C, we have 

N N 

2 = 1 2=1 

Therefore, we deduce that 

if, / M 1 GAr(Rei) 
sup \RNit)\ < sup I ... 

^ ,^supexp((^g[/,/] 

< (6-23) 
for some constant C"(To). By the two-constants lemma [NN22], (6.21)-(6.23) imply 

VT6]0,To[, sup |i?^(t)| ^ (C"(^o)w)-^^(^■^")/^ </>(r,ro) = arctan( ^ ° ). (6.24) 



In particular, for any compact K of the complex plane, we can find an open disk of radius Tq which 
contains K, and thus show (6.19) with i?Ar(t) e o(l) uniformly in K. □ 
We observe from the proof that Proposition 6.3 cannot be easily extended to Tjv \ t\ e 0(1) with 
Tat +00. Indeed, the ratio GAr(TAr(Ret))/|GAr(rArt)| in (6.23) wiU not be bounded when N ^ co, 
hence applying the two-constants lemma as above does not show i?jv(i) ^ 0. 

Corollary 6.4 In the model with fixed filling fractions e, assume the potential Vq satisfies Hypothe- 
ses 5.1. Then, if Lp : A ^ C can be decomposed as ip = ipi + p>2 with <p\,<pi extending to holo- 
morphic functions in a neighborhood of A, then the model with fixed filling fractions e and potential 
V = Vq + (p/N satisfies Hypotheses 5.1. Therefore, the result of Proposition 5.5 also holds: the 
correlators have a 1/N expansion. 

Proof. Hypothesis 5.1 contrains only the leading order of the potential, i.e. it holds for (Vo,e) iff it 
holds for {V = Vq + H/N, e). Proposition 6.3 implies a fortiori the existence of constants G+, G_, G > 
such that: 

G_ G^ s; I Z^l^^^ I s; G+ G^ (6.25) 
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Using this inequality as an input, we can repeat the proof given in Section 3 to check to obtain Corol- 
lary 3.7 (i.e. the a priori control reminded in (5.9)-(5.10)) for the potential V. Then, in the recursive 
analysis of the Schwingcr-Dyson equation of Section 5 for the model with fixed filling fractions, the 
fact that the potential is complex- valued does not matter, so we have proved the 1/N expansion of 
the correlators. □ 

This proves Theorem 1.3 in full generality. 

6.4 expansion of n-kernels 

We can apply Corollary 0.4 to study potentials of the form: 

K,x(a = ^ - E Hx, - (6.26) 

j 

where Xj e C\A, and thus derive the asymptotic expansion of the kernels in the complex plane, i.e. 
Corollary 1.7 and 1.8. Indeed, let us introduce the random variable Hc{x) = Z;"=i cj Xiili ^''^{^j ~ ^i)- 
We now know from Proposition 6.3 that In /x]^'^^[e*^'= is an entire function. Therefore, its Taylor 
series is convergent for any t e C, and we have: 

if„,c(x) = exp(lnA.];;^,3[e*«"=(-)]) 
1 r ^ AC 

= exp(2^f,n2S(2^^1"(^-'-^'))^-(^^'---'^'-0 (^-2^) 

r^l ■ i = l 7 = 1 



which can also be rewritten: 



where we introduced: 



i^„,e(x) = exp(2i/:f;[W^,]) (6.28) 



£c,x/(a;) = E / (6.29) 

j — l J 00 

As a consequence of Proposition -5.5, Wn e 0(iV^^") and has a 1/N expansion. Therefore, only a 
finite number of terms contribute to each order in the n-kernels, and we find: 



Proposition 6.5 Assume Hypothesis l.S. Then, for any K ^ —1, we have the asymptotic expansion: 

K , fc+2 



X„.e(x) = exp I f; 7V-^- ( ;^^x,c[Vt^i'=>]) I , (6.30) 



where 5 = inf j d(a;j , A) is assumed larger than 5n introduced in Definition 5.4- For a fixed K, it is 
uniform for x in any compact of (C\A)". □ 

If iy9 : A ^ C is a function such that 0{s) e o(|s|~'') when \s\ — > oo, one could study by similar methods 
the asymptotic expansion of the exponential statistics where H[Lp] = Xlili fi^i): that we 

would establish thanks to Corollary 6.2 up to o{N~-^^^^), where 

'-^%^^\ (6.31) 



2{k + 9) 

Note that K{r) > implies r> 1 + k + 29 = 7/2. In particular, we can deduce: 

Proposition 6.6 The central limit theorem 5.6 holds for test functions Lp such that |<^(s)| e 
Q(|g|-(i+K+2e)) ^^^^ ^ 
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7 Fixed filling fractions: 1/N expansion of the partition func- 
tion 



In this Section, we restrict ourselves to the fixed filhng fraction model, i.e. we study y^'^ ^ for some 
e e Eg. Following § G.2, it is again not difficult to consider potentials of the form F = Vi + V2, with 
Vi, V2 is holomorphic, so we will write down proofs only for holomorphic V . 

7.1 Interpolation principle 

Recall that, if (14)t is a smooth family of potentials so that dtVt is holomorphic in a neighborhood of 
A, we have: 



We are going to interpolate in two steps between the initial potential V, and a potential for which the 
partition function can be computed exactly by means of a Selberg /3 integral. 

7.1.1 Reference potentials 

We first describe a set of reference potentials. Let 7 = [7", 7"''] be a segment not reduced to a point, 
and p- two elements of {+1}. We introduce a probability measure supported on 7: 



da^.p(x) = ^Jix-j-)p-ij+ -x)P^ dx, (7.2) 

where the constant Cr,p ensures that the total mass is 1. It is well-known that a^^p = fjXq'"^ for the 
following data: 

• if (/9~, p~^) = (1,1), (Tj^p is a semi-circle law, and it is the equilibrium measure for the Gaussian 
potential 

7^ + 7+ \ 2 

(7+-7-)n-^ 2 J' (7+ -7-)^ 

on 7, any interval of 7™ = M which is a neighborhood of 7. 



if {p ,p^) = (—1,1), <7-y^p is a Marcenko-Pastur law, and it is the equilibrium measure for a 
linear potential: 

^ 4(.-7^) ^ 2 

"f'P ' "f'P V -^7 

on 7, any interval of = [7~, +oo[ which is a neighborhood of 7. 
if (p^, p+) = (1,-1), we have similarly a linear potential: 

4(7+ - X) 2 

on 7, any intevral of "f„i =] — co, 7+] which is a neighborhood of 7. 

if {p~ , p'^) = (—1,-1), (T-y.p is an arcsine law, and it is the equilibrium measure for a constant 
potential: 

V^,p = 0, c^,p = 1 (7.6) 

on 7 = 7 = 7™. 
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V '7 

When we choose 7 = 7™, the partition function Zj^'^p' of the initial model with such potentials are 
special cases of Selberg f3 integrals, and therefore can be computed exactly. We have in general: 

Z]:,y''- = exp ( - [m)N' + (1 - m)N] In (7.7) 
where -Z^^ do not depend on 7 and their expression is given in Appendix B.2. 
7.1.2 Step 1: interpolation with a reference potential 

Given A = U^^q A/j = U^=o['^h ' '^h ]' consider the support of the equilibrium measure mJ^''^, which 
is of the form Sg = |Jh=o ^h,e = [Jh=ol-^h e' '^h e}' '^^^h signs ph = {p^^Ph) hrdicating the soft of hard 
nature of the edges. Let {Uh)o^h^g be a family of pairwise distinct neighborhoods of Ah. We denote: 

V^ciA^) = J lu, (x) [en ^s....p. (x) + J] 2e„, / In |x - e|das„,^,p„ (O) • (7.8) 

h=0 h'lth •'Sfc/, 

By construction, K-ef.e is holomorphic in the neighborhood U = lJ^=o of A, and: 

a 

CTrcf,£(a;) = 2 e/icrs„.,,p^ (7.9) 

is the equilibrium measure for the potential K-of,e on A. Indeed, it satisfies the characterization (1.20). 
Notice that CTrof has same support as /Xeq,e, edges of the same nature, and same filling fractions. 
Besides, Kof satisfy the assumptions of 1.2. Then, if we consider the convex combination of potentials 
Vs = {\— s)V + sVrcf , it follows from the characterization (1.20) that the equilibrium measure associated 
to Vs on A with filling fraction e is precisely: 

& = (1 - s)pX^% + S fTref.c (7.10) 

Besides, since both pX(i% ''^'^^ CTref.e satisfy (5.2) with edges of the same nature, we conclude that if V 
satisfies 1.2, so does the family (K)s€[o,i]- Therefore, we can use Theorem 5.5 to deduce from (7.1) 
the asymptotic expansion: 



N,e,l3 



,.if--|f S/-'f^.ii(''-«)-KK))<«'^-«)|. (7.11) 

7.1.3 Step 2: localizing the supports 

We now have to analyze the partition function for the reference potential Vrcf defined by (7.8). When 
g = (the one-cut regime), K-ef coincides with one of the reference potentials, and we know that up 
to exponentially small correction, 2^'"^'^ will be given by the Selberg integrals described case by case 
in (B.6)-(B.8), so there is nothing more to do. 

Assume now g ^ 1. Let us define shortening fiows on the support: 

IJ , IT I Q <J 



if (Ph'Ph) = (I'l) (-I'-l)' we set S^^^ 

if (Ph^Ph) = (-1^ 1)' we set = [a^,e'"L + ^("/t,e - 

if (Ph'Pt) = (1' -1)' we set S^ ,, = [a+^ - <(a+^ - a^A,a+J. 



2 2 ' 2 ' 2 
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We consider the family of potentials on A: 

9 p 

Kla^) = ^ luja:) (e„Fs^,^,„Jx) - 2e,, / \n\x - ^\das^^ , (7.12) 

for which the equihbrium measure is obviously 2j^=o ^'i ''^s^ .ph ^'^'^ support S* = [Jk-o^h e- 
Accordingly, T^^f ^ satisfies Hypothesis 5.1, uniformly for t in any compact of ]0, 1]. Besides, the 

partition function 2^^°^'^ can be computed exactly in the limit t — > 0. If we introduce: 

f K/,. + <J/2 if {pl,pt) = {lA)ov{-l,-l) 

<H={ ^Ih if {PH,Pt) = {-^^) , (7.13) 

I <H if iPH,Pt) = (^^-^) 

we find that: 

— n i<.-<.'r^""'> (7.14) 

where is the maximum allowed interval associated to which depends on the nature of the edges 
(i.e. on Ph) as described in § 7.1.1. We remark that the dependence in t factors and: 

zl'l"/^'' = exp f^-[(l3/2)N'el + (1 - /3/2)iVe.] In (^^^^^^) | Z'n%,,^, (7.15) 

where Z^^^ is an analytic function of Nch- The asymptotic expansion when ^ of those factors 
associated to reference potentials is described in Appendix B.2. We just mention that it is of the form: 



Zjj'f^ = Ar(^/^)^+^^ exp ( 2 , (7.16) 

fc>-2 

for 7 = [7^,7^]. Therefore, we obtain: 

^kS = n 1""'^ - "° '^'i' fl ^n\.,p \-[wm'{ t el) + (1 - mN_ 



In 



^^^ + (1-^/2)^ 
s s 



X exp I 2 N-'^ f'dsUm)(t'iy 

+ / ^ (c'.K^t J(0 [ • (7-17) 



,k^-2 h=0 

Is, 

By construction, the integrand in the right-hand side is finite when s — > 0. The expression does 
not make it obvious for the terms fc = — 2 and fc = — 1, but it can be checked explicitly since the 
eigenvalues in different S| ^ decouple in the limit s — > 0, in the sense that: 

<"''(^) =, £ <'^""'(^) + 0(1), (7.18) 

and the expressions for the non decaying contributions to Wi when N is large are given in Ap- 
pendix B.l. 
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7.2 Expansion of the partition function 

We establish in Lemma B.l that the partition functions for the reference potentials Z'^ ^ do have an 
asymptotic expansion of the form: 

= Ar(/J/2)Af+ep exp ( ^ iV-'^ + 0(7V-=°)) , (7.19) 

where: 

3 + /3/2 + 2//3 (3/2 + 2//3 -1 + ^/2 + 2//? 

7++ = ^ : e+_ = e_+ = , e__ = ^ . (7.20) 

Therefore, we have proved a part of Theorem 1.3: 

Proposition 7.1 If {V,e) satisfy Hypothesis 1.1 and 1.3 on A (instead of B), we have: 

Zn'Xp = ^(^/')^+^ exp ( 2 N-"" F^^^ + 0(7V-=^)) . (7.21) 

with a universal constant e = Xjfi=o ^Ph depending on (5 and the nature of the edges. 

Let e, the equilibrium filling fractions in the initial model m]v'^- order to finish the proof of 
Theorem 1.3, it remains to show that the stronger Hypotheses 1.2-1.3 for /i^'^ imply Hypothesis 5.1 
for the model /i]^''^^ ^ for the model with fixed filling fractions e e £g close enough to e« , that all 

( — 21 

coefficients of the expansion are smooth functions of e, and that the Hessian of Fl with respect 
to filling fractions is negative definite. This last part is justified in Proposition A. 3 proved in Ap- 
pendix A.l, whereas to prove the first part, we rely on the basic result also proved in Lemma A.l and 
Corollary A. 2 in Appendix A.l: 

Lemma 7.2 If V satisfies Hypotheses 1.2-1.4, then (V, e) satisfies Hypotheses 5.1 for e e Eg close 
enough to e*. Besides, the soft edges a* and W"'g^^(x) are functions of e, while the hard edges 
remain unchanged, at least for e close enough to e* . 

We observe that, once wl~^^'' and the edges of the support a* are known, the W^} for any n > 1 
and A: > are determined recursively by (5.51)-(5.49) and (5.70)-(5.73), where the linear operator 
/C~^ is given explicitly in (5.21)-(5.27), and thus depend also analytically on e close enough to e*. 
Similarly, F^^^ for fc 5= are obtained from (7.11)-(7.17), which shows their analyticity for e close 
enough to e«. 

Corollary 7.3 IfV satisfies Hypotheses 1.2-1.4, then WnJ and F^ areC^ functions of e e Eg close 
enough to e^. □ 

This concludes the proofs of Theorem 1.3 and Corollary 1.7 announced in Section 1.4. 



8 Asymptotic expansion in the initial model in the multi-cut 
regime 

8.1 The partition function 

We come back to the initial model A*]^'^, and we assume Hypotheses 1.2-1.4 with number of cuts 
(g + 1) > 2. We remind the notation N = {Nh)i^h^g for the number of eigenvalues in Ah, and the 
number of eigenvalue in Ag is A'o = iV — I]fi=i ^h- The Nh are here random variables, which take the 
value Ne with probability ^^'^ p/ Z^'^. We denote €^, the vector of equilibrium filling fractions, and 
N^, = Ne^. Let us summarize four essential points: 
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We have established in Theorem 1.4 an expansion for the partition function with fixed filhng 
fractions: 

= iV(^/^)^+^ exp ( 2 ^^iS) ' (8-1) 



iV! 



where e are independent of the filling fractions 



• By concentration of measures, we have established in Corollary 3.8 the existence of constants 
C, C" > such that, for N large enough, 

A.]^;^[|Ar-Ar,| >lniV] sce^^'"^-*^'^'"'^. (8.2) 

• Thanks to the strong offcriticality assumption, we have after Lemma 7.2 that ^ is smooth 
when e is in the vicinity of e, . From there we deduce that, for any K,k ^ —2, there exists a 
constant Ck,K > such that: 



K-k (i^i^*)^-''^ 
7! 



^Ck,KN-^^+^'^\N-N,f-''+\ (8.3) 



• We establish in Proposition A. 3 in Appendix A.l that the Hessian {F^ „ ')" is negative definite. 



We now proceed with the proof of Theorem 1.5. 

8.1.1 Taylor expansion around the equilibrium filling fraction 



By the estimate (8.2), we can write: 



|Ar-JV»|«lnAf 

with: 

r^<(5+l)^e-S^^'"^^. (8.5) 

And, we have, for any K ^ —2: 

m 



0!iNi,...,N^siN 
\N.-N,\^\nN 

K K-k (pi'^hU) 

E exp ( £ 2 A^-('=+-') ./ ■ {N - iV,)®^' + iV-(^+i)i?K 

\N-Nt\^\nN 

And, since A^^(^+i)^J^ 1 for N large enough: 

|e^"'""''«- - 1| < 2\N-^^+^^Rk\ (8.7) 
^ N-(^+^) 2 2Cfe,K(lniV)^-^-+i <C^iV-(^+i)(lniV)^+3. 

fe = -2 

where we finally used (8.3). Notice that, since e^, is the equilibrium filling fraction, we have (F^^^^)'^ = 
0, and therefore, for any K ^ 0: 

K_ K_~k (^{fe})^) 
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exp ( 2 2_ iV-('^+^-) ■ {N AT.)®^") ( 

K 

1 + 2 A^-'=Ti5[Ar- AT,] + 0(Ar-(^+i)(lniV)^+3)), 



„ iTT T-, , 3 ■ ( AT - JV, ) ® ^ + 2 i 7ri> , , 3 ■ ( iV - AT, 



fc=l 
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where we have introduced: 



2i7r 



2m 



and for any vector X with g components: 



mi,...,m^>— 2 



(8.9) 



(8.10) 



Since the number of lattice points N satisfying |7V — 7Vj,| ^ IniV is a 0(ln A^), we can write: 



K 



Nl,...,NgEZ' 



fc=l 



+0(iV-(^+i)(lniV)^+4). 
where we have set 77 jv = In TV. 



.11) 



8.1.2 Waiving the constraint on the sum over filling fractions 

Now, we would like to extend the sum over the whole lattice Z^. Let us denote A*_^ 
mm Sp {-Fl ^ ')" > 0. For any a > small enough, there exists a constant C" > so that: 



JVeZ'' 
|JV-JV.|>>7jv 



|JV-JV,|>r)jv 

C" Volg(n)(n+l)V,-^*'^(i-")f"', 



(8.12) 



where Volg(n) = (2n+l)f-(2n-l)9 s: g29 n^'^ is the number of points in so that n |Ar-Ar^| < 
n + 1. Therefore: 



^ gi7r-r,,^-(JV-JV»)®^+2i7r^.,3-(JV-JV.)^jY - AT*)®^' 



-l+J g-->i*(l-a)s'ii;]V 



(8.13) 



where C^j is a constant depending on j. In other words, by unrestricting the sum in (8.11), we only 
make an error of order 0(e-^*(i"^)'), which is 0{N-^^). Then, we remark that: 



We have thus proved: 
rzVA ( K 







{v^,p\t.,p). (8.14) 



2i7r-' 







(-u,,^|t,,^) + 0(7V-(^+i)(ln7V)^+*). (8.15) 



The term appearing as a prefactor of A^^*^ is bounded when iV — > cx). So, by pushing the expansion 
one step further, the error 0(iV~(^+^)(ln A^)^+^) can be replaced by 0(iV~(-^+^)). This concludes 
the proof of Theorem 1.5. 
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8.2 Deviations of filling fractions from their mean value 

We now describe the fluctuations of tlie number of eigenvalues in each segment. Let P = {Pq, . . . ,Pg) 
be a vector of integers such that P — Nt.i,^h e o{N^I^) when N oo. The joint probability for h e [0, g\ 
to find Ph eigenvalues in the segment A/i is: 

f^NA^ = P] = TT9 5-T r^v-A (8-16) 



We remind that the coefficients of the large N expansion of the numerator are smooth functions of 
P/N. Therefore, we can perform a Taylor expansion in P/N close to e*, and we find that provided 
P — Ne^, e o{N^^^), only the quadratic term of the Taylor expansion remains when N is large: 







—1 ^ 

(v^M-r^.p)) exp [( |] (/3/2)(P;, - Ne.,h)') Iun] 



X exp (l (Fj-^V • (P - Ne.)®^ + (i^j-i*)' • (P - Ne.) + o{l)) (8.17) 



In other words, the random vector A7V = (AA^i, . . . , ANg) defined by: 

g 



AN^ = N,- Ne^j, + (8.18) 

h' = l 



converges in law to a random discrete Gaussian vector, with covariance 

[(fJ/V]"^- We observe 
that, when /3 = 2, fj^^^ = so that N — Ne^ converges to a centered discrete Gaussian vector. 

A Elementary properties of the equilibrium measure with 
fixed filling fractions 

A.l Smooth dependence 

In this section, we prove Lemma 7.2, i.e. we establish under some assumptions that the equilibrium 
measure /z^ ^ depends smoothly on the potential V and the filling fractions e. Let V be at least 
and confining on A. We know that the support consists of a finite union of pairwise disjoint segments: 



g 



Upon squeezing A, we can always assume that it is the disjoint union of (g + 1) pairwise disjoint 
segments A = [jfi^Q^^h, which are neighborhoods of S^^ in M. The Stieltjes transform wj ^^'^ of 
this equilibrium measure satisfies: 

^helO,gl ^xeSl^, + iO) + + iO) = y'(a;), (A.2) 

and 

V/i6lO,5l, S ^ = e/.. (A.3) 

The general solution of (A.2) takes the form: 
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where PY is a polynomial of degree g which should be determined by (A. 3), and we recall the notation: 

This implies that the equilibrium measure has a density, which can be written in the form: 

dfi^ix) = dxSYix) fl\x al,r'^. (A.6) 

for some p* = + 1, and SY{x) is a smooth fimction in a neighborhood of the support. This rewriting 
assumes 5^ (a*) 7^ when p* = —1. 

Let I = 10,5] X { + !}, and V. = {I e I pi = —1}. We assume strong off-criticality as defined in 
Hypothesis f.2: 

Hypothesis A.l The initial data V\Q\ and e[0] is strongly off-critical, in the sense that, for any 
I e I\n, 5(a/) ^ and S'{ai) ^ 0. 

For any (h, •) e TL, we consider the values aj fixed and equal to a'^^j'"^^- We introduce the open set: 
U = |a 6 Y\ ^-ll ct/ 1 foi" e I are pairwise distinctj, (A. 7) 

and the map: 



defined by: 



(F,r) : X Rg[X] X C^(A) X RS+l , ^2g + 2-\Hl ^ ^g+l 



Vlen Fi[a,P,V,e\ = , + U) ?r ■. 

' ^ 2m aj — X 



V/iG[0,gl n[a,P,V,e] = -e^ + (j) ^w[a,P,Vm, (A.9) 



A, 2i^ 



where: 



u;[a,P, KJ(a:;j = -| r -,/ ^ + ^ ^ TfTT ' (A.IO) 

a[a](x) = flVa; a/. (A.fl) 

lel 

By construction, the data of the equilibrium measure jj^ satisfies: 

{F,T)[a^,PY,V,e] = 0. (A.f2) 
We would like to apply the implicit function theorem to show that: 

Lemma A.l If Hypothesis A.l holds and V[0] is C (resp. analytic) with r > 2, cx^ and PY are 
Qr-i ^j.ggp^ analytic) functions of {e,V) close enough to (e[0\,V[0\) . 

The contour integrals (^^ used in this Appendix A.f can be rewritten as integrals over the segment 
S. Thus, these manipulations do not assume that V is analytic in a neighborhood of A, and it does 
makes sense to consider only V of class C in the statement of Lemma A.f. 
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Proof. To achieve this we need to show that {da,dp){F,T) when evaluated at a point such that 
(F, T)(a, P, V, e) = 0. We first compute: 

It is well-known that the g x g submatrix of (A. 13) consisting of indices h,h' e ll,gj is invertible 
when a e U, and on top of that, we find: 



^h'elO,gl Ji^[a,P,y,e] = = _ Res 

^ cA'* /a 2i7r cr[Q;J(^) «^gc 



e dC 



^. (A.14) 



In particular, this expression does not vanish for h' = g, hence (A. 13) is invertible. Then, we compute: 
[a, P,V] = Sjjf^ ""t"]^^^ - «/) - (V'iai) - V'iO) 



da J ' ' ' Ja 2i7r (a/ - ^) 

F/[a,P,F]-Fj[a,P,l/] 



+ 



2{ai - aj) 

P\aj) a[«](OdC V"{ai){£,-ai)-{V'{ai)-V\Oy 



Sij(- 

= lim V^^(54?«[a,P,F](a;)), (A.15) 

and the latter does not vanish when evaluated at (a^^Qp pY^^\ ^[0], e[0]) thanks to Hypothesis A.l. 
So, (da,dp)(F,T) is invertible at this point. Besides, if V is C for r > 2, then d(F,r) is C''^. 
Therefore, by the implicit function theorem, for (e, V) close enough to (e[0], 1^[0]), there is a unique 
function function (a^, P^) so that {F, T){a.Y , P^ , V, e) = 0. By uniqueness, it must correspond 
to the data of nYq,e- D 

Corollary A. 2 If Hypothesis 1.2 holds for (y[0],e[0]) with V\Q\ not necessarily real-analytic hut 
with r 5= 2, they hold also for (V, e) close enough to {V\Q\, e) and V\Q\ C, and the density (A.G), once 
multiplied by (j{x), is for such data. 

Indeed, (A. 4) shows that W^^^^'^ (x) is a linear function of Pe^^\x)/a{x) and hence, once multiplied 
by a{x), it is a smooth function of e. □ 

A. 2 Hessian of the value of the energy functional 

We are now in position to prove: 

Proposition A. 3 If Hypothesis 1.2 holds for (V^[0],e[0]) with V[0\ not necessarily real-analytic but 
at least , pj ^''^ is for (e, V) close enough to (e[0], V^[0]), and the g y< g matrix with purely 
imaginary entries: 

\fh,h' 6 [l,gl, {tY)h,h' = ^ (A.16) 

2i7r OEhOCh' 

is such that ImrJ^ > 0. 

Proof. We are going to justify that Fe is a concave function of e in the domain: 

{ee]0,ip+\ 2e/. = l}, (A.17) 

h=0 
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and within Hypothesis 1.2, which will imply the result. Let e a vector of such that 2^=0 = 
and e ^ 0. For any e e £g, we define e* = e + te, which belongs to £g for t small enough, and 
/J,* = /J,^ . Its characterization as an equilibrium measure imply that : 

U\x) = V{x) - 2 y In |x - y\dfi\y) (A.18) 

is locally constant on the support of /i*, and /u*(A/i) = e/i + tch- Let us denote C/^ the value of U*{x) 
when X e Sh- 

Let us denote Ff,(t) = F^t After a classical result (see e.g. [AG97, Joh98]): 

F.it) =^[jj\^\x- y\dii\xW{y) - j V{xW{x)\ (A.19) 

It follows from Corollary A. 2 that the density of /j,* is away from the edges, and its derivative is 
integrable at the edges, thus define a measure that we denote i/*. It has same support as /i*, and 
we deduce from the characterization of that x i— > —2 Jlnja; — y\dv^{y) is locally constant on the 
support of and i^*[A/i] = e/j. We deduce that Ff,{t) is and: 

K{t) = ^ / (2/ \r.\x-yW\y)-V{x))dv\x) = -^j]uUh (A.20) 

•' h=0 

eh does not depend on t, whereas 1 1-^ Ul is as one can deduce from the expression (A.18). Thus, 
Fe is and: 

F':{t) = /3 2 / In \xh - y\d,y*{y) en (A.21) 
where xu is any point in the interior of A/j. This can be rewritten: 

F'm = pj] [[ 2\n\x-y\d,y'{x)l;,,{x)diy\y)U,{y)=f3f]Q[iy'U„iy'Uj (A.22) 

h^Q-'-' h=0 

It is well-known property (see e.g. [AG97, Dci99]) that, for any signed measure v with total mass 0, 
i/] 5= 0, with equality iff = 0. Since we chose e ¥= 0, the vector of measures (z^*lA,Jos;h$g is 
not identically zero, hence F^{t) < 0. In other words, F^ is strictly concave for any direction e, hence 

1—21 

F^ ' is a strictly concave function of e e £g so that Hypothesis 1.2 holds. □ 

B Model Selberg integrals 

B.l Non decaying terms in correlators 

Let us denote Wi''''^ the first point correlator in the model with potential Vp^^^^ described in § 7.1.1. 
Let us denote A = (7+ - 7")/4. It admits a 1/iV expansion: Wip,jm = Y,k^-i ^''^ W^^'^''^"' '' .The 
expression for the equilibrium measure gives access to Wi ^^■'•''"'''^ 

T^{-i};7™,++. N _ a;- V(x-7-)(x-7+) 



Ix- 




X — 




Ix- 


7"' 



= _L(i_^/^_^), (B.l) 



Wry^^-^ix) ^ (B.2) 

wt'^-'^-~^{x) = ' (B.3) 

V(a;-7 )(a;-7+) 
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and then we deduce from (5.49): 



Besides, we find from (5.51) that ah these models share the same W. 



{0}. 
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W^-^-ix^x,) = I ^ f - 1 + x,x.-(x, + x,)ll±^+7-7^ ^ ^ 
/5 2(xi-X2)2V V(^i - 7-)(a:i - 7+)(^2 - l-){x2 - 7+)^ 

B.2 Exact formulas for partition function 

Let us denote: 



. «p{42f!.(:-f)f]4f)},,,-n™. ,B.e, 

r ^ 
= n |A.-A,|^n^"'^^-^^''^dA„ F_+(A) = A 

l^i<j^Ar i=l 

2^:^ = / n iA,:-A,fndA^ 

^ (r(i + (j-i)/3/2))'r(i + j/?/2) 

y r(2 + (iv-2 + j)/3/2)r(i + /3/2)- ^ •''^ 

These are the values of the reference partition functions given in (B.6)-(B.8). To emphasize that they 
can be defined for N not restricted to be an integer by analytic continuation, we introduce a function 
related to the Barnes double Gamma function: 

r2iN;h,h) = exp(^ C2{s;h,b2,x)), (B.9) 
Vds s=o / 

where: 

1 e~*^ t*~^ dt 

Us;h,b2,x)=^J^ (1_,-M)(i_e-M)- (B.IO) 

Its properties are reviewed in [Spr09], in particular it solves the functional equation: 

r2(x + 62; 61,62) = j?^|^V2^6f-"/'\ r2(l; 61,62) = 1. (B.ll) 

We deduce from (B.ll) the representation: 

n r(l + J/3/2) = (27r)^/^(/3/2)/^^V4+^(i/2+m) ^^J(^^+l)^^^^ (B.12) 
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Therefore, we can recast the Selberg integrals as: 

Z++ = exp ( - il3/4)N^ lnN+ {l3/4 - 1/2)NItiN + 7V[(/3/2) ln(/3/2) + ln(27r) - lnr(l + ^/2)]) 

r(jV+l) 
r2(iV + 1 ; 2/(3, 1) 

Z+-^ = e:ip - {l3/2)N^lnN + (13/2 - l)NliiN + N[pin{/3/2) + Iti{2tt) -liiT{l + 13/2)]^ 

r^(iV+ 1) 

''r{l + N(3/2)Tl{N+l; 2/(3,1) 

Z^-p = exp(^/3iV2ij^2 + Ar[(2-/3)ln2 + (3/3/2)ln(/3/2) + ln(27r)-lnr(l + /3/2)]) 

T{2/(3 + N - 1)T{N - 1) V^{N + 1)T2{2N - 1 ; 2//3, 1) 



r(2//3 + 2N - l)T{2N - 1) r2(l + NP/2)Tl{N + 1 ; 2/(3, l)T2{N - 1 ; 2//3, 1) 

B.3 Large asymptotics of the partition function 

We need the asymptotic expansion of Barnes double Gamma function [Spr09]: 

irr 0IR^^ ipx^ l( B^ 3+^/2 + 2/^ 
lnr2(x;2//3,l) =2:^00 \ ^ ^ 2 2 / ~ 
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-X'(0; 61,62)+ ;2(fc-l)! 62) x-'= + 0(x-=^), (B.13) 

Ek{bi, 62) are the polynomials in two variables appearing as coefficients in the expansion: 

1 



(1 - e-''i*)(l - e-''2t) 



=^ 2 Ek{b^,b2)t'' (B.14) 



fc>-2 

which are expressible in terms of Bernoulli numbers. x(s ; 61, 62) is the analytic continuation to the 
complex plane of the series defined for Res > 2: 

X(5; 61,62)= ^ ( k 1 h ^s - (B-15) 

For instance: 

x'(0;l,l) = -^ + C'(-l) (B.16) 
We remind also the Stirling formula for the asymptotic expansion of the Gamma function: 

where Bk are the Bernoulli numbers. We deduce the asymptotic expansions: 
Lemma B.l 

lnZ++ = -{3(3/8)N'^ + (/3/2)iV In TV + ( - 1/2 - (3/A + {(3/2) ln(/3/2) + ln(27r) - lnr(l + (3/2))N 

+ ii^!^lniV+V(0:mi) + !^f:l + o,l, (B.18) 

In Z+-^ = -(3/3/4)iV2 + (/3/2)iV In iV + ( - 1 + (/3/2) ln(/3/2) + ln(27r) - In r(l + /3/2)) iV 

8/2 + 2/(3 ,^ , fii(/3/2) ln(27r) , , , , 

+ ' g +2x'{0;2/(3,l)-^^ + ^^ + oil) (B.19) 

InZ^^ = (/3/2)iVln7V + ( - /3/2 + (/3/2) ln(/3/2) + (/3/2 - 1) ln(2) + ln(27r) - lnr(l + /3/2))iV 

+ IniV + 3x'(0 ; 2//3, 1) + - 13(2/g - 13(/3/2) ^^^^^ 

-ln(/3/2)+ ^^+0(1) (B.20) 
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where the o(l) have an asymptotic expansion in powers ofl/N. 
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